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In this thesis massive higher derivative gravity theories are analyzed in some detail. One- 
particle scattering amplitude between two covariantly conserved sources mediated by a gravi- 
ton exchange is found at tree-level in D dimensional (Anti)-de Sitter and flat spacetimes for 
the most general quadratic curvature theory augmented with the Pauli-Fierz mass term. 
From the amplitude expression, the Newtonian potential energies are calculated for various 
cases. Also, from this amplitude and the propagator structure, a three dimensional unitary 
theory is identified. In the second part of the thesis, the found three dimensional unitary 
theory is studied in more detail from a canonical point of view. The general higher order ac- 
tion is written in terms of gauge-invariant functions both in flat and de Sitter backgrounds. 
The analysis is extended by adding static sources, spinning masses and the gravitational 
04 Chern-Simons term separately to the theory in the case of flat spacetime. For all cases the 

microscopic spectrum and the masses are found. In the discussion of curved spacetime, the 
masses are found in the relativistic and non-relativistic limits. In the Appendix, some useful 
calculations that are frequently used in the bulk of the thesis are given. 
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Part I 



Introduction 



I. 



INTRODUCTION 



General Relativity (GR) [T] can be thought as a modification (albeit a major one) of Newtonian 
gravity. The history of gravity starts with Newton's ^ force law in 1687. The precise form of this 
law is 



where G is the universal Newton's constant and rm are the masses of the interacting particles. 
The minus sign indicates that the force is attractive. The Celestial mechanics gave successful 
verifications of this law. A well-known example of the success of Newton's theory is the prediction 
of a new planet. After the discovery of the orbit of Uranus there were attempts to calculate 
this orbit theoretically. However, these attempts did not work exactly since there were deviations 
between the calculated and observed orbits [2j [3j. Using Newton's law ([T]) the location of Neptune 
was predicted and just after this calculation it was observed exactly at that location [21 E]- Hence 
an "outer" mass was responsible for the perturbations of the orbit of Uranus. Newton's law, 
with the help of "dark matter" (Neptune in this case) worked well for the case of Uranus. After 
this resolution, there were studies of all planets. When Mercury was studied, another problem 
was found: The calculations of precession of Mercury's perihelion was in contradiction with the 
observed value. The calculations gave larger value than the observed ones. The first attempt to 
solve this problem was again to invoke the "dark matter" idea. It was thought that there must be 
a planet, Vulcan, in the solar system j3lll]. However, this did not solve the problem since Vulcan 
was never found. The solution came with modifying gravity: replacing Newtonian gravity with the 
Einsteinian one [3]. 

In GR the spacetime geometry is determined by the matter, energy, pressure etc. and the 
dynamics of matter or light is determined by the geometry of spacetime. GR is constructed 
with the idea of equivalence principle and the general coordinate invariance (and with the added 
assumption that the equations are wave- type equations with second derivatives on the basic fields) . 
The Einstein equation is 



where the left hand side is the geometry part and it is known as the Einstein tensor and the right 
hand side is the matter part. In (|2]), R^^y and R are the Ricci tensor and Ricci scalar, respectively. 
These quantities can be thought as a measure of the curvature of spacetime. They depend on 
the spacetime metric g^i^ and its derivatives. The connection is taken to be metric compatible 
(V^ = 0). In the right hand side, there is a coupling constant related to the Newton's constant 
which in four dimensions is k = SvrG, and T^,/ is the energy-momentum tensor coming from the 
matter sector. 

Einstein's equation, right after was introduced, was immediately applied to solve the then more- 
than fifty year old problem of Mercury's orbit. After the Schwarzschild solution was found, which is 
the unique static spherically symmetric time-independent solution to ([2| , the perihelion precession 
of Mercury was calculated and the answer was found to be compatible with the observed value 
13 EHZ] • There were also other predictions that came out from GR such as the bending of light 
passing by the sun and the gravitational redshift of light which were also supported by observations. 



F = - 



Gmim2 



(1) 



(2) 
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Inspite of these successful solutions and predictions there are some observations that GR in its basic 
form as given in ^ cannot explain. The data taken from supernova explosions fSHTU] show that 
the universe has an accelerated expansion. GR cannot explain this phenomenon in its pure form^. 

If GR is taken as the correct theory, there must exist a dark energy component in the matter- 
energy budget of the universe which can be represented as a constant term, the cosmological 
constant, A. The density of the dark energy is estimated to be 75% of the total energy density of 
the universe and numerically it reads 

PA = 4x10-6^, (3) 

which is very small compared to the vacuum energy density that comes from the energy density of 
the vacuum of quantum field theory (QFT); 

PA ~ 10~'^^pQFT, (4) 

|1L)H12| . This mismatch of the "experimental" and "theoretical" values of the energy density of the 
vacuum is known as the cosmological constant problem. Therefore, it can be thought that GR is 
not the full theory and it must be modified in the infrared (IR), namely at ultra-large distances or 
ultra- weak interacting regimes, such that without a need of dark energy the accelerating expansion 
of the universe happens. [Admittedly, this alone will not solve the question of why QFT vacuum 
is almost empty.] One of the suitable candidates for this modification is to give a tiny mass to the 
graviton in the theory, in such a way that massive GR still passes the solar system tests. 

From the perspective of QFT, at low energies, GR can be thought as a weakly interacting 
massless spin-2 field. The story starts with the classification of particles with respect to their spins 
and masses, as degrees of freedom (TUl US]- The representations of these degrees of freedom are 
fields. For spin-0 mode of massless particles the representations is a scalar field, and for spin-1 
it is a vector field. For spin-2 mode the symmetry is the general coordinate invariance when the 
interactions are taken into account, and around flat and maximally symmetric spaces, one has a 
massless helicity-2 particle. This points to GR [lUI [TSH^O] with the action 

I=-[ (fxy/^R , (5) 



where g is the determinant of the spacetime metric and the rest are defined in section V B There- 
fore, it is also natural to combine the quantum theory with GR. 

There are other motivations to try to modify GR and perhaps construct a quantum gravity 
theory. The Schwarzschild solution of (|2]) is a black hole with a curvature singularity, to explain 
or resolve this type of singularities a quantum gravity theory is needed. Also, to understand the 
very early universe when both the quantum effects and gravitational effects are dominant at the 
same time, a unification of gravity and quantum mechanics is necessary. There are attempts to 
construct quantum gravity theories. String theory is one such candidate. Higher curvature gravity 
theories also are candidates. With this motivation it is also possible to add higher curvature 
terms to Einstein gravity ([s]). The higher curvature terms are negligible at low energies, but they 
dominate at high energy domains. Another reason to introduce the higher curvature terms is that 
the Einstein-Hilbert action ([5| is not renormalizable. Higher curvature terms make the theory 
renormalizable but ruin the unitarity, namely yields negative norm states in the scattering matrix 
1211 122]. 



^ Tfj,i, is taken as the standard energy-momentum tensor of the matter. 
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To summarize the argument, we can say that GR needs modification at both the IR and 
ultraviolet (UV) regimes. Therefore, in this thesis Einstein gravity will be modified both by adding 
mass and higher curvature terms. Hence we modify the theory at both the UV and IR regimes. 

The outline of the thesis is as follows: In the next sections of this chapter the massive gravity 
and higher derivative gravity models are discussed separately. In the massive gravity model mainly 
the Pauli-Fierz [23| [M] model is studied and also the nature of the van Dam-Veltman-Zakharov 
[251 ESj discontinuity between the strictly massless theory and the arbitrarily small mass theory 
is discussed. Then some specific three dimensional massive theories are studied. After that, 
some technical details will be given about the relevant spacetimes, the cosmological constant and 
the higher dimensional gravity models. At the end of this Chapter the higher derivative Pais- 
Uhlenbeck oscillators are briefly reviewed. In the second chapter the general higher curvature 
massive gravity theory is analyzed. Its unitarity structure is studied by computing the tree level 
scattering amplitude in generic D dimensions. Also, the Newtonian potentials are calculated. 
The third chapter is devoted to the analysis of the three dimensional unitary theory. These two 
chapters are based on the papers ^^Massive Higher Derivative Gravity in D- dimensional Anti- 
de Sitter Spacetimes" j27] and Canonical Structure of Higher Derivative Gravity in 3D " |2H] 
respectively. Then, the conclusion part comes. Also, an appendix part is added so that some of 
the calculations can be followed easily. 



II. MODIFYING GRAVITY 

In this part the modification of the Einstein theory is discussed. By modifying gravity we 
mean that we change the degrees of freedom of the theory in some way. First the mass term is 
added to the theory which changes the degrees of freedom from two to five in 3 + 1 dimensions. 
Also, adding higher derivative terms change the degrees of freedom. The massive gravity theory 
is discussed first. Then the higher derivative gravity theory is considered. At the end of this part, 
two 3-dimensional massive higher derivative gravity theories are introduced. 



A. Modifying Gravity with Mass Terms: 

There are two ways to give mass to a gravity theory. One is to add directly the mass to 
the action. The other way is to introduce scalar fields. When these fields are evaluated at the 
background the general covariance is broken and the graviton gets mass. Both ways end up in the 
same class of theories. In this discussion we will follow the first way. The mass is added to the 
Einstein-Hilbert theory in a Lorentz invariant way as follows 



I 



-g 



1 



-R 



1 



(6) 



which is known as the Pauli-Fierz (PF) action, the second part is the PF mass term |23| [M] and 



Tj^i, is the 
- /i^y, and 



m" is the mass parameter. We first consider this theory in a flat background. In ([6 
flat spacetime metric and h^,y is the linear part of the metric perturbation, g^y = rj^^ 
h = rj^yh^^ is trace of this metric perturbation. To see more explicitly how the added PF part 
gives mass to the graviton, we should linearize the theory around the background spacetime. From 
both Q and ^ , the linearized Einstein tensor can be written in the absence of sources as 



0, 



(7) 
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where = d^d^ and the hnear parts of the Ricci tensor and Ricci scalar are 

If the metric perturbation is constrained to be transverse and traceless, that are d^h^y = and 
/i = 0, then yields 

a^ (V - r]^.^.h) = 0, (9) 

and for massive particles it is known that the Klein-Gordon equation must hold that is 

- m^) </> = 0, (10) 

for the mostly plus signature of the metric. To write (|9]) as a massive field equation, the PF term 
must be added: 

{d^ - m^) (V - Vi^uh) = 0. (11) 

Since this term cannot be obtained from the curvature terms, it is added to the action ([s]) by 
brute force. Later on, it will be seen that modifying GR with higher curvature terms will generate 
this mass term automatically when the coupling constants of the higher curvature terms have 
a special combination. From the transverse and traceless conditions it can be seen easily that 
the PF theory describes a massive spin-2 field. The metric perturbation is a symmetric rank-2 
tensor, so that it has ten independent degrees of freedom in 3 + 1 dimensions. However, from the 
transversality condition there are four constraints which eliminate four degrees of freedom leaving 
six. The tracelessness condition also eliminates one degree of freedom and the theory has five 
degrees of freedom. We must also note that this mass term is the unique ghost and tachyon-free 
Lorentz-invariant combination. When the sign is changed in the middle of the PF mass term, then 
it produces tachyons in its propagator structure [TD] . 

The massless limit of the PF theory must tend to Einstein gravity as the usual continuity 
arguments in physics dictate. However, the rm? — )• limit does not give the results of the m? = 
theory. For two point sources the interaction potential, in the Newtonian limit of GR, in four 
dimensions is 



However, the massive theory gives 



C/ = -^^. (12) 



t/ = -'^^, (13) 
3 r 



where the Newtonian potential is greater than the usual one. This effect is known as the van 
Dam-Veltman-Zakharov (vDVZ) discontinuity^ |25l I26j . The linearized PF theory is different from 
the linearized GR in the massless limit. This discontinuity appears in flat backgrounds |29H31] . 
However, in curved spacetime this discontinuity does not appear. From the general amplitude 



equation (|79|) it can be seen more explicitly. There is a ^ fraction by which the flat spacetime 



limit and massless limit do not commute. Going to flat spacetime limit first the vDVZ discontinuity 
appears in the massless limit, but taking the massless limit first the discontinuity disappears in 



Redefinition of the Newton's constant does not solve the problem, since then the defiection of light changes by 
25%. 
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four dimensions |27| I3U] . Hence, the introduction of a cosmological constant to the theory solves 
this problem since the smallness of the mass can be compared with another measurable quantity. 
[One could argue that discontinuity has been mainly replaced by the non-commutativity of the 
limits, which is a valid objection.] From the pole structure of the general amplitude equation (79 1, 
it can be seen that for de Sitter spacetime (A > 0) a pole produces a tachyon and it is absent for 
anti-de Sitter (A < 0) spacetime^. 



B. Modifying Gravity with Higher Curvature Terms: 

To make GR perturbatively renormalizable, quadratic curvature terms aR^ + PR^^R^^ + 
^RficFvpR^'^'^'' must be added to the Einstein-Hilbert action (|5| [21]. But adding these terms 
makes the theory nonunitary because of the non-decoupling ghost introduced by the middle term 
j3Rfj_uR^''" [22^ . When this term is excluded the unitarity of the theory is regained but the theory 
becomes nonrenormalizable. To have a consistent quantum gravity, both of these properties must 
be provided. 

Before discussing further, the notion of unitarity must be explained. With unitarity, it is meant 
that the theory must be both ghost and tachyon free. Ghost is a particle that has a negative kinetic 
energy and tachyon is a particle that has negative mass-squared (in flat space) . In a more technical 
language the signs of the propagators must be correct that is (for the (— , -|-, -|-, -|-) signature) the 
(scalar propagator) 

D (p) ~ , ^ , (14) 

where p is the four-momentum. In the canonical form, the Lagrangian should be 

£ = ^^(t,x)(n-m2)v(t,x), (15) 

the kinetic term is the d'Alembertian operator and it must be positive. The mass-squared must 
again be greater than zero. Note that only the tree-level unitarity is considered here. The Feynman 
diagram of the tree-level scattering amplitude is shown in the figure . 



C. Topologically Massive Gravity 



In the discussion to follow, we will need a specific 3D theory that has been studied a lot in the 
literature. This is the topologically massive gravity (TMG) theory introduced in j33| [51] with the 
action 



I 



-R 



2/i 



3 



(16) 



Here ^ is the coupling constant of the gravitational Chern-Simons term. The sign of ^ is not fixed. 
^x^ll' jg ^YiQ three dimensional anti-symmetric tensor which is connected to the Levi-Civita symbol 
as e^^^ = — ;^/=e'*'^'^- The pure Einstein-Hilbert action does not propagate any degrees of freedom 
in 3D. After the gravitational Chern-Simons term is introduced, TMG propagates a single massive 



^ We should also note that Vainshtein claimed that [321 the discontinuity disappears at the nonlinear level once the 
finite Schwarzschild radius of one of the scattering particles is taken into account. 
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Figure 1: The tree- level one-particle scattering amplitude between two background covariantly conserved 
sources mediated by a graviton exchange. 



particle with helicity +2 or —2 (not both since the theory is not parity- invariant). The action for 
this theory can be written in terms of the gauge-invariant functions^ (134) as 



1=1 I d^x 



- ((t>q + a'^) + -a{q+a(f)) 



(17) 



Taking the variation with respect to the q field one of the other fields can be eliminated. This 
variation yields 



1 1 

-(j) + -a = 0, 

K H 



(18) 



from which 



-^(T follows. Putting this result back into (17) yields the action 



n 2 



(19) 



From (19) it can be seen that this theory propagates only a single massive degree of freedom with 
mass m? = Not to have a negative kinetic energy term, k must be chosen negative which yields 
the "wrong-sign" Einstein-Hilbert term. In fact this model propagates a single massive scalar spin- 
2 mode [S3]. Since the spin of the particle depends on the sign of fi, this model is a parity violating 
theory, as noted above. The detailed discussion and calculations for this model is given in Chapter 
3. 



■* Note that the gauge-invariant functions will be discussed later in the third chapter. 
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D. New Massive Gravity and Critical Gravity 

Another interesting higher derivative massive gravity theory is the recently introduced "New 
Massive Gravity" [33 [30]. The action for this model is 



If... ^ r „ ^ 

m 



d^x^g\-R+ — {Rl^- R^]\, (20) 



where there is a relation between the coupling constants of the higher curvature terms as 8a+3/3 = 
and the mass reads in terms of the coupling constants as m? = For the unitarity of the theory, 

K must be negative, and /? > 0. This model also implies the wrong-sign Einstein-Hilbert term. 



This model is not a higher derivative theory at the linearized level. The canonical form of (20) 



is (see more details in section XIII ) 



"M-^r^+s^i^j+ir'+s^'Jl' '''' 

where it can be seen that this model has two modes propagating with the same masses. Therefore, 
this theory is parity invariant. Also, in this form the necessity of negative k can be seen explicitly. 

After this theory was introduced [2S] a research activity started in massive 3D gravity theories. 
Especially in Anti-de Sitter (AdS) backgrounds some classical solutions and conserved charges 
of this theory was given j27| I5^46| . The question of how the theory gains mass was answered 
in |47| US]. In this work Higgs mechanism of the general quadratic theory in three dimensions 
was given by writing the theory in a Weyl invariant form. Breaking the Weyl invariance in the 
background, graviton gains mass, scalar field remains massless and the Weyl gauge field gains mass 
or can remain massless depending on the parameters. 

The four and higher dimensional extensions of this theory, called "the critical gravity", was 
also investigated in |49H51] . For these cases the special points are again set to zero 3a + /3 = 
and 4a{D — 1) + Df3 = for four and the generic D dimensional cases respectively, for which the 
massive scalar mode vanishes |27j . In the critical gravity, the mass of spin-2 is set to zero and 
the energy of black holes become zero for the AdS vacuum. It is also found that the energy of 
the massless spin-2 modes vanishes on-shell. There remains only spin-2 logarithmic modes with 
positive energy. 



III. CONSTANT CURVATURE SPACETIMES 

In our discussion we will need the maximally symmetric backgrounds which are the flat, Anti-de 
Sitter (AdS) and de Sitter (dS) geometries. Maximally symmetric means that the curvature of the 
spacetime is constant everywhere and in any direction. Therefore, if the curvature is known in one 
point it is enough to know the curvature at every point. There are only three maximally symmetric 
spacetimes that are Minkowski, which is flat so that the curvature is zero everywhere, dS and AdS. 
These spacetimes are classified with respect to the curvature scalar R, the dimensionality of space- 
time and the signature of the metric. De Sitter space is a {D — l)-dimensional surface embedded 
into a D dimensional flat spacetime for which the signature of the metric reads (— , -|-, -|-, -|-). 
For the AdS spacetime the signature of the embedding metric reads (—,-|-,-|-, ...,—). For AdS 
spacetime there is no horizon but for dS there is a horizon [521154] . Examples of the maximally 
symmetric flat spaces are planes with Euclidean signature (-I-, -|-, -|-) and its higher dimensional 
generahzations. Spheres (^S^^ are dS or positively curved spacetimes, and hyperboloids are 
negatively curved or AdS spacetimes. 
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The Riemann tensor of a maximally symmetric D-dimensional spacetime is 

R 

(22) 

where the curvature scalar R is constant and D is the dimensionality of the spacetime [7]. For 
AdS or dS the curvature scalar is R = (fji^- Here, A is the cosmological constant and for A < 
the spacetime becomes AdS, negatively curved spacetime, and for A > it is dS, positively curved 
spacetime. If A = then the spacetime is flat. 

The cosmological constant is related to the vacuum energy density in D = 4 as 

Pvac = (23) 

Here p is the energy density of the vacuum which is proportional to the pressure of a perfect 
fluid Tfj^u = (p + p) UfJJu + pgpu- To see this proportionality the energy- momentum tensor in the 
Einstein equation is divided into matter and vacuum pieces. 



Rpu - ^g^^uR = SttG (t^, + r;r , (24) 



1 

2' 

where the vacuum energy- momentum tensor is T™"^ = —pva.c9^iu- The vacuum energy-momentum 
tensor becomes a perfect if — Pvac = Pvac- The Einstein equation can be written with a cosmological 
constant term as 

Rpu - \gpvR + Ag,.u = SnGT^u- (25) 

The relation between the vacuum energy density and the cosmological constant can be seen from 
these equations [7j . The maximally symmetric solutions of this equation are dS and AdS (and also 
Minkowski) spacetimes, for the pure vacuum case. 

In the first part of the analysis the higher curvature massive gravity is studied in a D dimensional 
(A)dS spacetime (AdS/j). In this way the effects of higher dimensions can be seen. Also, with this 
general approach a unitary theory may be found in D > 5. 



IV. CANONICAL ANALYSIS 



In this part some technical details of the canonical analysis are given. To put a Lagrangian 
into the canonical form means to write it in terms of scalar fields that are in the form of harmonic 
oscillator Lagrangians. The field equations of these Lagrangians are Klein-Gordon equation. So 
that, KG equation can be seen as simple harmonic oscillators. From these oscillators the particle 
spectrum can be studied. Also, for investigation of the unitarity of a theory, this form is useful. 
Therefore, first the KG equation for free particle is reviewed. Then the generalization to higher 
derivative of this equation known as Pais-Uhlenbeck oscillator is discussed. 



A. The Klein-Gordon Equation 



The story starts from the quantization of the non-relativistic energy equation for a free particle. 
The energy of a free particle is 



(26) 
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where E is the energy and p is the momentum of the particle. To quantize this, the variables are 
promoted to being operators by taking E ^ H = ihdt and p —ihV, where h is the reduced 
Planck constant, and ^ = dt = do, V = ^ for one dimensional space^. With these substitutions 
(26) turns to the well known Schrodinger equation 

ihdtip {t, x) - 



2m 



(27) 



where ip (t, x) is a complex function. A relativistic version of this equation can be obtained by 
upgrading the relativistic dispersion relation^ |55) 



E^ 



2 2 I 2 4 

p c + m c , 



to an operator equation 



?/.(t,x) =0, 



(28) 



(29) 



where c is the speed of light. Going to natural units c = 1, h = 1 and defining the d'Alembertian 
operator □ = 9^9^, where the components of the partial derivative is d'^ = (^d^, (in Minkowski 
background with the signature (— , +, +, +)) yields'^ 



0, 



(30) 



which is the classical relativistic wave equation of a massive particle. Taking the Fourier transform 
of ip{t,x) = J ^^e^^'^ip {t,p) and putting it into (29), the equation of motion for a harmonic 
oscillator is found with the frequency ujp = \/\p\'^ + m^, that is [551 EZ! 



ij{t,p) = 0. 



(31) 



In the field theoretical approach the same result can be found by taking the real spin-0 particle 
Lagrangian density with mass m 



C{t,x) 



(32) 



where the field depends both on time and space ip = il){t,x). The Euler-Lagrange equations 
- daTrM-Tx of (I32I) gives (Isol [Ml EZl • 



of (32 



9V '"I'didfi;) 

After getting some basic informations about the relativistic wave equation, the higher derivative 
terms can be added to the Lagrangian to generalize this discussion. The field equations of the higher 
derivative Lagrangians are known as the Pais-Uhlenbeck oscillators | i58j . 



B. Pais-Uhlenbeck Oscillators 

The Lagrangian density of the higher derivative real scalar field is |59j 

^=-^V'(n(°+"^'))v'. (33) 



^ Here the derivative is df_i = (J^, V) and the four momentum is P** = (E, p) = id'^ in the natural units. 

^ If the equation is not squared the operators will be inside the square root E = \J p^c? + vnP-d^ — > ihdttp = 

\J (— iftV)^ + m'^c'^ip. This form has some problems in evaluating the operators. Also, it is nonlocal. 

Note that for the signature (+, — , — , — ) the KG equation reads (□ + m^) ip = 0. 
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For simplicity N is set to two, i. e. = 2, and the examination is done in the nonrelativistic hmit, 
that is all the space derivatives are dropped. For these conditions the Lagrangian density becomes 



(34) 



where cji is used instead of rrij. Also, note that the wi = W2 case will be different from the 
tui 7^ u}2 case, where the first case is the degenerate case and the Lagrangian density becomes 
purely quadratic 



C = - 



The Euler-Lagrange equation for ( 34 ) is 



(i2 /dC\ d fdC\ dC 



dt"^ \ dq ) dt \ dq ) ^ dq 



(35) 



(36) 



which can be obtained by taking the variation of (34) with respect to the variable q = q{t). For 
n-th order variational problem ( 36 ) takes the form 

/ dC 



N 
n=0 



(37) 



in the non-relativistic limit. This equation (37) can be obtained by taking the variation of the 



general Lagrangian L = L (^q, q, q, . . . , q^"-^^ , which is [BDl EI] 
5S 



/ dtSq 
Jt, 



N 

■E 

L n=0 



dy dL 
Jt) (9g(") 



+ 



N-l 

E 

Ln=0 



E ^.("X 



(38) 



where Pq(n) is defined as 



N 

E 

k=n+l 



d 
dt 



k—n—l 



dL 



(39) 



The Hamiltonian for the generic case is written fiSCTl lUT] as 



N-l 



n=0 



From ( 37 ) the equations of motion can be found for = 2 as 

g(^) + (uf + ujl)q + ujfujlq = 0, 



(40) 



(41) 



where q^'^^ is the fourth time derivative of q. The Hamiltonian (40) and the momentum (39) of the 
system for N = 2 reads 



H = -U<f-2q'^^^-{u^l + 



OJr 



2\ a2 



2 2 2l 

■uj^U}2q j 



and 



Pg(o) = (w? + u;^)q + q'^^^ and P^w = -q, 



(42) 



(43) 
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respectively. For the non-degenerate case coi ^ uj2 with the Pais-Uhlenbeck variables Qi = q + 
and Q2 = 9 + can be diagonalized as 

where the second term is a ghost term. However, for the degenerate case there can be a unitary 
region for proper constants [59] . 

V. CONVENTIONS 

In what follows we will use the following conventions: 

A. Flat Spacetime: 

The flat spacetime metric has the mostly positive signature d\a,g{ri^u) = (— )+,+;+,••• ) in 
its diagonal components. Since the spacetime is flat, the curvature tensors and scalars are zero 
R'^p^ = R^^ = R = 0. The partial derivative reads = {d\ -d') = ((9°, -V) and = 

{—dt, —di) = c?o, ~^)- 'Th^ d'Alembertian operator becomes = 3^ = d^d^ = —Oq + df = 
— (9o + V^. The Greek indices run as // = 0, 1, 2, . . . that counts all the coordinates of the spacetime. 
The Latin indices run as i = 1,2, 3,... and they denote the spatial coordinates. 

B. The Curved Spacetime: 

The metric has the signature (—,+, +,+,...) and the Christoffel connection reads 

= la"^ idf^9ux + dyg^,x - dxg^,u) , (45) 
and the commutation relation of the covariant derivatives are 

[V^,V,]V'^ = R^^\V\ (46) 

where 

^tifx = 9fj,T1x - duTJ^x + ^Jip^tx - ^Ip^'^^x (47) 
is Riemann tensor. The covariant derivative acts on the covariant and contravariant vectors as 

v^uj, = d^iv,-T^^,Lox, v^F^ = a^y^ + r;:;,y\ (48) 

The Ricci tensor and Ricci scalar are defined as 

Rf,u = g"" Rpf.au R = g^''Rf,u, (49) 

respectively. The d'Alembertian is □ = V^V^. The Greek and Latin indices have the same 
meaning as in the flat spacetime case. 
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Part II 

Massive Higher Derivative Gravity in 
D-dimensional Anti-de Sitter Spacetimes 

VI. INTRODUCTION 

In gravity there is no theory which is unitary and renormahzable at the same time. To get 
a renormahzable theory in four dimensions, higher derivative terms, aR^ + f^Rf^^, are added to 
Einstein-Hilbert action [21] . In doing so one can gain renormahzabihty but loses unitary. The cou- 
pling constant of the square of the Ricci tensor introduces a non-decoupling ghost term. However, 
if we omit this term then we gain unitary yet lose renormahzabihty. If a theory is non-unitary 
it shows itself as a repulsive force in the Newtonian limit between static sources. Therefore the 
theory has a better UV behaviour because of this repulsive force. In field theory this is what usu- 
ally happens: To have a better behaved theory, ghosts are introduced during the renormalization 
process but they decouple at the end if the theory is unitary. Therefore, bartering unitarity with 
renormahzabihty can not be accepted. 

In three dimensions there is a perturbatively renormahzable and tree-level unitary theory |35| I36| 
in flat background for a special choice of coupling constants of higher derivative terms that is 
8a + 3/3 = and for a reversed sign of Einstein-Hilbert term |62l - [55] and this theory is named 
as "New Massive Gravity" (NMG), that is a parity-preserving spin-2 theory. However, we do 
not know that if this special ratio between a and /3 will survive renormalization at a given loop 
level. Another interesting thing about this theory is that at the linearized level the theory has 
a massive graviton with helicities ±2 in its spectrum. Therefore, the theory has a non-linear 
extension to the Pauli-Fierz (PF) mass term. For this reason this theory solves an old problem 
of finding a non-linear extension of PF mass term. The equivalence of NMG and PF was shown 
in [SS]®. The physical meaning of this equivalence must be well understood since NMG theory 
is background diffeomorphism invariant not only at non-linear level but also at the linear level. 
However, PF theory is invariant under Killing symmetries of the 2-1-1 dimensional Minkowski 
spacetime. To have a better understanding about these symmetries a quite interesting approach 
was put forward in [BS] : In the absence of the Einstein-Hilbert term, the Weyl invariant form of 
the linearized NMG is written. Therefore, introducing the Einstein-Hilbert term at the linearized 
level breaks this symmetry and produces the mass of the graviton. As a result one can think that 
Einstein-Hilbert term provides the mass and the higher derivative terms gives the kinetic energy. 
This perspective explains the sign change of Einstein-Hilbert term. Looking back, this result is 
expected since pure Einstein gravity is non-dynamical and does not give any propagation in three 
dimensions. It is like the mass term in a scalar field theory at the linearized level. When a kinetic 
energy is introduced to the theory then it plays a role in the dynamics. In the case of NMG the 
kinetic energy comes from the higher derivative terms. Consequently, the mass of the graviton 
comes from Einstein-Hilbert term by breaking Weyl invariance. This point may be important for 
constructing massive gravity theories in other dimensions. 

In this chapter, the most general quadratic curvature gravity theory will be considered in a D 
dimensional (anti)-de Sitter background. We also augment this theory with a PF mass term. We 
will study the propagator structure of this theory by finding its one-particle exchange amplitude 
between two covariantly conserved sources. For this reason we first linearize the action 



Note that, beside this parity-preserving theory, there is also the parity violating Topologically Massive Gravity 
(TMG) IMIEI]. 
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1 = 1 d''xV^!^^R-'^ + aR^ + /3Rlu+7{Rluap-^Rlu + R')] 

+ 1 d«xV=5|-^(/i^.-/iV^matter}> (50) 

where Aq is the bare cosmological constant, k is related to the D-dimensional Newton's constant 
as K = 2^£)-2Gd where is the D — 2 dimensional solid angle. The other parameters are the 

coupling constants a, /3, 7 and is the mass parameter. In total we have a seven parameter theory 
which is the most general quadratic model including the dimensions and cosmological constants^. 
Because of these parameters the theory potentially has various interesting limits and discontinuities. 
At this level, there are no constraints on the parameters. They may be positive, negative or zero. 
After computing the tree-level amplitude we will constrain the parameters not to have any ghosts 
or tachyons in the theory. Also, some terms do not contribute to the equations of motion for some 
specific dimensions. For example, for D = 4 the Gauss-Bonnet term becomes a surface term and 
for D = 3 it vanishes identically. Therefore, in three dimensions all the information is carried by 
the Ricci tensor and the Riemann tensor has no more information. In three dimensions the model 
can be extended by adding the Chern-Simons term /i (tOT -\- |r^^ to (50). But here we will stick 

to (50). In the next chapter we study this case in detail^*^. The theory reduces to model with 
the PF term for two dimensions. In this chapter we consider D >3. The theory that we consider 
has general covariance except for the PF mass term. 



Various limits of the spin-2 model which is defined by the linearization of (50 ) have been studied 
in the literature. Though, for some certain limits of the above action there may still be some 
interesting new models. One of them is NMG which is the case for D = 3 and in flat background 
without the PF mass term and 8a + 3/3 = 0. Apart from finding such new models we will also 
explore discontinuities of the full seven parameter theory. The discontinuities come out while the 
order of limits are changed when some of the parameters approach zero. One such discontinuity 
is the so called van vDVZ discontinuity. The resolution of this discontinuity comes from the 
introduction of the cosmological constant. Then taking ^ — >■ limit GR results are recovered at 
tree- level |i29Vl32| [HH]^^. However, the discontinuity reappears when the quantum corrections are 
taken into account [7D]. Up to now the linear theory has been considered. Once, we consider the 
non-linear theory a ghost arises in massive gravity and this effect is known as the Boulware-Deser 
instability |71j . 

The ingredients of this chapter is as follows: In the second section we will analyze the linear 
equations of motion which is obtained around an AdS background. We also discuss certain special 
limits. In the next section we calculate the one-particle scattering amplitude and the following 
sections will be devoted to the discussions of some limits and discontinuities. At the last section 
the conclusions and discussions will be given. 



In the absence of the source terms and at the hnearized level, one can reduce the number of parameters in the 
action [661 167) , but here for the sake of generality we shall work with ( 50 1. 
See [591 168) for this case, without the higher curvature terms. 

Another resolution of the discontinuity may follow even in flat space if the Schwarzschild radius of the scattering 
objects is taken as a second mass scale in the theory |32) . 
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VII. LINEARIZED EQUATIONS 



The field equations can be found by taking the variation of the action (50) with respect to the 
metric. After variation we get the field equations as follows: 



1 



+27 



R 



1 



1 



g^yR + Kog^y + 2aR R^^ - -g^,uR + (2a + (3) {gf^^n - V^V^.) R 



+PU R^ 



1 \ f 1 \ _ 

-gpuR ) + 2/3 ( Rfiaup - ^gpuRap ] R"'' + — {hpu - gpuh) = t^u- (51) 



Here r^i, is energy-momentum tensor coming from the source terms, etc... . The background 
metric (jf^i, namely the vacuum, is a non-singular solution to the field equations in the absence of 
the matter terms. The vacuum is the (anti)-de Sitter space which is the maximally symmetric 
vacuum with the Riemann, Ricci tensors and Ricci scalar curvature given respectively as 



R 



2A 



ppucr 



^9pv9pa 9p.a9vp) 1 R) 



2A 



:9iJ.u, R 



2DA 



(52) 



We made all the contractions with respect to the background metric. If we put the background 



metric (52 ) in the field equations (51 ), the effective cosmological constant A can be found in terms 



of a, (3, 7, and Aq: 



A-Aq 
2k 



+ /(a,/3,7)A2 



(53) 



where /(a,/3,7) = {aD + (3) + l ^^^xj^_2) • This is a second order equation with respect 

to A. Therefore it has two solutions which are 



A 



4k/ 



1 ± + 8k/Ao 



(54) 



and for reality of the vacuum 8k Aq/ > —1. Hence the spacetime can be both dS or AdS. One 
of these solutions vanishes in the absence of Aq. In this case the non- vanishing root becomes 
A = — 2re7' There are some exceptional points of this equation of motion for which the effective 
cosmological constant becomes equal to the bare cosmological constant. In four dimensions A = Aq 
since /(q,/3,7) becomes zero. Also, in three dimensions if we set 3a + l3 = then we again 
have f{a,f3,^) = and A = Aq. Moreover, if we set 7 = and aD + (3 = 0, which gives 
/(a, /3,7) = 0, the theory again has Aq = A. In three dimensions without setting 3a + f3 = 0, we 



have Ah 



l±y/l-8K{3a+(3)Ao 
4k(3o+/3) 



and for A to be real 1 > 8k (3a -|- (3) Aq. 



The next step is to find the linear equation of motion of (51) around the constant curvature 
background, ^^jy = g^^ + h^^. By using (54) we eliminate the bare cosmological constant and after 



linearization the equation of motion becomes [72] 



Tf,, {h) 



2A 



g^, + (2a + /3) ( g^,n - V^V, + ^^g^.u ] R 



+ 13 (ug^, - ^—f9puR^\ + 1^ (V - 9puh) 



M2 



(55) 



19 



Note that the PF term is already hnear. Here we define a new constant in terms of the other 
coupling constants as 



1 AKD 4A ^ 4A(D-3)(D-4) 

+ {D-l){D-2)^- 



(56) 



T^u {h) is the energy-momentum tensor which contains all higher order terms and the source 



Tay. Q!^j^ is the linear cosmological Einstein tensor 



2A , 
D-2 



(57) 



and are the linearized Ricci tensor and the linearized scalar curvature R^ = {g^'^R^ 



respectively. They read 



Oh 



VaVuh] , R 



2A 
D-2 



h. (58) 



To calculate the scattering amplitude we need the trace of (55 1, that is T = g^yT^^^ and it reads 



UaiD - I) + DI3)U - ID ~ 2) [ - + 4/A 



R' 



M2 



[D - l)h = 2T. 



(59) 



This equation is a wave equation, and for Aa{D — 1) + Dfi = something special happens, since 
the dynamical part of this equation vanishes. Actually, for D = 2> this choice will lead us to the 
NMG. After computing the tree level scattering amplitude between two sources, and constraining 
the theory to have no ghosts or tachyons will give us this special point. Before moving on we can 
still analyze some limits of the theory at the linearized level. For this purpose we drop the sources. 
Since there are no sources we can not get the unitarity regions but we can capture the ranges of 
the parameters in which tachyons drop out. The discussion bifurcates whether is zero or not. 
First we will consider ^ 0. 

A. Massive case: 



We take the divergence and double divergence of (55) to find the constraints on the deviation 



part of the metric, hny Since, (55) is divergence free without the mass term, the constraints come 



from the divergence of the mass term as 



vv^ V - n/i = 0. 



(60) 



Replacing the covariant derivatives that appear in the background Ricci scalar in ( 58 ) which is 



2A 



{D-2) 



2A 



{D-2) 
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h. For the flat space case we have = 0. 



and using the first equation in ( 60 ) lead to = — 
This choice forces the field to He traceless (59), and from the first equation of (60) it becomes 
transverse. The field /i^i, has ^^^^^ independent components. The transverse-traceless condition 
gives D + 1 equations that eliminates + 1 of the independent components. Thus we are left 
with the remaining {D + 1)(D — 2)/2 independent components. The linear Einstein tensor becomes 



d hp^ui where we have changed the places of the covariant derivatives in (58). With 



these constraints the field equation for the remaining independent components of the field becomes 



0, 



(61) 



which is a quadratic equation that describes two massive excitations with masses 

1 1 



2k/3 2|k/3| 



1 + 4K/3M2. 



(62) 



To have real masses > —4^, and at the saturation point the two masses become equal. These 
masses become non-tachyonic when the parameters are chosen properly. If k and /3 have the same 
sign, one of the excitation becomes negative and produces a tachyon. If they have opposite signs, 
the term inside the square root must be < 4 | k/3 | < 1. In this case both excitations become 
non-tachyonic. However, this theory is not unitary as we will see. 
For the generic A case the trace equation reads 



(4a(L> - 1) + Dl3)n -{D-2) 



/i = 0. 



(63) 



In this equation we can see that h becomes a dynamical scalar field. The dynamical part can be 
eliminated by setting Aa{D — 1) + Dfi = 0. After choosing this special point we still have two 
options to satisfy (63). We may set /i = and the field becomes again transverse and traceless. 



and the field equations read 



(a + /?□) 



2A 



[D-l){D 



2 



+ 



M2 

2k 



0, 



(64) 



and it has generically two excitations 



2A 



{D-l){D-2) 2/3 



2A a y M2 

XD-l){D-2)^2p) 



(65) 



After computing the one-particle exchange amplitude we will again discuss this theory and see that 
for special points, this theory will give us a ghost and tachyon free theory. The other option to 



satisfy (63) is tuning the mass as 



D-1\k 



D - 1 



{D-l){D-2) 



(66) 



which is the partially massless point, vanishes for fiat background and for curved background a 
higher derivative gauge invariance appears [75] . Therefore, the field has one less degree of freedom 
(DOF) compared to the massive one. Moreover, the mass is allowed to be negative in AdS as long 
as it satisfies the Breitenlohner-Freedman type bound. In four dimensions higher derivative terms 
vanish and the mass only depends on the cosmological constant, = In higher dimensions 
the partially massless theory depends on the higher derivative terms. 
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B. Massless case: 



Without the Pauh-Fierz mass term the theory becomes invariant under background diffeomor- 
phisms 5(^h^y = V^^j^+Vi,^^, since = and 5^R^ = 0. As we mentioned above the divergence 

and double divergence of (55) do not give any constraint on h^j,u. Therefore, becomes a dy- 
namical variable for T = 0. If the coefficient in front of the d'Alembertian operator is fixed to zero 
than the dynamics of vanishes. The equation (59) is satisfied for two conditions: R^ can be 



zero or the coefficient, ^ + 4A/, can be set to zero in which case R^ need not vanish. The later 
one is not acceptable because the field equations leave a gauge invariant object undetermined. 



VIII. TREE-LEVEL AMPLITUDE 



To find the Newton potentials and to analyze the particle spectrum of ( 50 ) , we are going to find 



the tree-level scattering amplitude between two covariantly conserved sources which is defined as 



A 



1 



d^x ^--gT' {x) h'^'^ix), 



(67) 



by considering the full theory ( 55 ) . In this equation T^^^ is one of the sources and h^^'^ is the 
deviation which is produced by the unprimed source. The factor ^ is put to get the correct 
Newton's constant. Since the field has independent components we need to decompose h^'^ to 
eliminate the unphysical parts of it. Therefore, the decomposition must be done in such a way 
that the physical components of /i^;^ will be determined by T^'^. The usual choice is to define 



(68) 



where /i^J is the transverse and traceless part of the deviation. is the vector part with a 



symmetrization that is defined with a ^ factor and it is divergence free, (j) ^iiid V are scalar 



functions. With this definition (68), the amplitude equation becomes 



A 



1 



d'^x 



-g(T'^X^^^ + T'i; 



(69) 



where the terms in the middle of ( 68 ) vanish since we have covariantly conserved sources and the 



total derivative vanishes at the boundary by use of Stoke's theorem. In (69 ) the tensorial quantities 



and Tp must be written in terms of the trace of energy momentum tensor. In order to get this we 



take the trace, divergence and double divergence of (68) 



h = \Il(j) + Dtlj, 



2A 



{D-2) 



0(1) + n^p, 



(70) 



where we used VV^/i^i^ = Oh. This condition comes from the nonzero mass term and it is not 
a gauge condition. From these two equations we can write 01;^ in terms of Oip. We hit the first 
equation with the □ operator and subtract it from the second one which gives us 



{D -1){D- 2] 
2A 



□V, 



(71) 



which can be put into the first equation in ( 70 ) yielding 



(72) 
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This equation can be inserted into ( 59 ) so that ip can be written in terms of the trace of the energy 
momentum tensor as follows 



= J - + 4A/ - cAD 



M2 
2k 



(D-l) 



[D-l)iD-2)^ 



2A 



n + D\ T, 



(73) 



where c = ^^^_^2°' + 25?^- Now we will write /i^J such that it is determined by the trace of the 
energy momentum tensor. First of all, we will find the transverse and traceless part of the field 



(2) 

equations by using the Lichnerowicz operator acting on spin-2 symmetric tensors 

A^L^hf^u = —Hhfj^i, — 2R^puah^" + 2R^ (p,^y)p- 
Some properties of this operator (that we need) were collected in 



(74) 



,(1) 



(2), 



Using these properties we have 



A 



(0). 



2A 



{D-2 



-h 



TT 



(75) 



(76) 



With this equation the transverse traceless part of the deviation can be written in terms of the 
transverse traceless part of the energy momentum tensor as 



'TT 

''/lU 



2 (/3n + a)(A^^ 



(2) 



D 



4A , 

) + 



" fiu ■ 



(77) 



The only thing that is left is to decompose the energy momentum tensor so that we can connect 
h^J with the energy momentum tensor. To find this we first decompose the energy momentum 
tensor as (68). Then we take the double divergence of this equation, keeping in mind that we have 
covariantly conserved sources. Also taking the trace, one can find 



D 



-T + 

1 D-l 



V V + '^-^ 



X □ + 



2Ai:» 



{D-l){D-2) 



-1 



T. 



(78) 



We are ready to compute the scattering amplitude by putting (73 77 [78 ) into ( 69 ) . After doing 
manipulations by using ( 75 ) and suppressing the integral sign the amplitude equation becomes 



AA = 2r;,<j(/3n + a)(Ai^ 
2 



(2) 



4A 

D-2' ~V 



rpflU 



+ 



D - 1 



T' {{pn + a){n + 



4A . 
D-2' 



+ 



4A 

{D-2){D-ir 
2 



r {{i3n + a){n + 



{D-2){D-1) 



T' {- + 4A/ - cD 



4A 
D-2 

M2 



(79) 



) 



M2 
K 



-1 



2kA 



(D-l) 



□ + 

□ + ■ 



2AD 

:D-2){D-1) 

2AD 

D-2){D-l), 



T. 
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From this main result we can consider various limits. Solving this integral for non zero cosmological 
constant is highly nontrivial. However, the particle spectrum of the theory can be considered by 
analyzing the pole structure of the amplitude. For the general case there are four poles which read 
as 

- 2AD 

□i = 

□2,3 = -A-r^ + T^,]^\ir^ + T^,] -/3(77F^- — )^ (80) 




□4 = 

To have tachyon free theory these poles must be positive. By choosing the constants properly 
such a theory can be found. However, this is not the only restriction on the parameters that 
appear in the general quadratic theory. With these poles the residues can be found. These residues 
must be negative in order not to have any ghost terms. This is the other restriction that we use 
while finding a ghost and tachyon free model. However, in the most general form the residues are 
cumbersome. Therefore, we will restrict the theory and compute the poles and residues for some 
interesting limits. Also, we will compute the Newtonian potentials for these limits. 



IX. MASSIVE THEORY IN FLAT SPACETIME: 



Looking at (79), we can see that the massless limit and the flat space limit do not commute 
because of the 



-j^ term. Therefore, if we want to get the usual Newtonian potentials in flat space 
and without mass term, the flat space limit must be taken before the massless limit. Otherwise, 
we encounter the well known vDVZ discontinuity. Let us first look at the limits which produce the 
vDVZ discontinuity that means we first take the flat space limit A — )■ and go to massless limit 

M2 ^ 



0. With the flat space limit (79) becomes; 

-1 



AA = -2T'Apd^ + -() 



D - 1 



T' \ i3d^ + -d' 



M2 



(81) 



(2) 

here the Lichnerowicz operator goes to = —3'* and this can be seen from (74). The spectrum 



of (81) has two massive excitations which can be found by solving the second order equation for 
d'^ and the masses are the same as (62). Modifying (81) with the mass terms give 



I3{mi 



T' 



1 



1 



Q2 



rrij 



Q2 



rpflU 
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{D-1 



1 



1 



Q2 



nij 



Q2 



m_ 



T 



(82) 



For P < and /? > 0, (82) produces a massive ghost. Therefore, /3 must be set to zero to avoid 



ghost terms. The Newtonian potential energy (U) can be calculated by use of Green's function 
technique. For simplicity we take the sources as Tqq = mi6(x — xi), Tqo = m2S(x' — X2), where mi 
and 7712 are masses of the sources, and the other terms of the energy-momentum tensor are taken 
to be zero. Moreover, we calculate these potentials in four and three dimensions. The Green's 
function for modified Helmholtz equation — k"^ reads ^-^"0(^1^1 ~ i"2\) and '^^||,^'_^^| |74| for 
two and three space dimensions, respectively. Then the potential energies read 



1 



m\'m2 



u 



2fS{m'^-m'i) 47r 
mir?i2 ]_ r. 
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D = 3, 
Z) = 4, 



(83) 
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where r = \xi — a?2 1 and are defined in ( 62 ) . If we take the /? — >■ limit to take care of the 



ghost term, (83) becomes 



U 



4Gmirn2 —Mr 



D = 3, 
D = A. 



(84) 



The first equation of (84) is the Newtonian limit of massive gravity in three dimensions and it 
gives attractive force when k is positive. For this case the M — )• limit does not exist. As a; — >■ 0, 
Kq{x) — >■ — ln(x/2) + 7£;, which gives the expected ^ force for small separation of the sources. 



The second equation of (84) is the Newtonian limit of massive gravity in four dimensions. Unlike 
the three dimensional case, M — )■ limit exists. However, the | term indicates the famous vDVZ 
discontinuity. G is the Newton constant and it is taken as Therefore, in three dimensions 

massive gravity gives the correct Newtonian limit despite that in four dimensions it does not give 
the expected limit. 



X. MASSLESS THEORY IN FLAT SPACETIME: 

In the above discussion we see that first taking the flat space limit and then going to massless 
limit does not give us the expected Newtonian limit. In this section we first take the massless limit 



and then go to the flat space limit. When we set = and A — in (79 ), we get 



AA = -2r;,|/3a^ + ia2| r^- + ^^^^r'|/39^ + ^s^j r 

2 -T'{cd^--dA~\. (85) 



{D-l){D-2) I n 



Unlike (81), (85) has three poles that are 



S? = 0, 81 = -2^. 8| = 1^ (86) 

To make a full analysis we need the residues of these poles which read 

i?e.f52) _ 2^ 2k(£-2) / ,x _ 2^ 

From the second pole we see that k/3 < for not having tachyon and form the residue of this pole 
K < not to have a ghost. For negative k the residue of the massless pole fixes the dimension to 
three, since for > 3 makes the residue positive in which case it produces a massless ghost. The 
residue of the third pole becomes positive for negative k. To eliminate this residue one has to set 
c = 0. For these conditions the theory becomes unitary and has its special name as New Massive 
Gravity (NMG). If we look to the Newtonian potential which is written as 

U = ■^mim2 {Koinigr) - Ko{mor)) D = 3, (87) 

87r 

where = — ^ and uiq = ^^ga+sp) • negative k, the second term gives a massive ghost which 
decouples in the case of 8a + 3/3 = 0. In this case only the first term lives and gives an attractive 
force. Also, if we choose the Pauli-Fierz mass term as M = mg the NMG has the same Newtonian 
limit as the usual massive gravity. For D > 3, a massive ghost comes out and does not decouple 
unless ^ = 0. For example we can look at D = 4: 

U = (l - ^e-™«^ + ) , (88) 
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where = 2k{^+i3) • '^^^ term in the middle is the ghost term [2P. When we take /? — )• the 
ghost term vanishes and taking nia — >■ oo the last term drops out and we are left with the usual 
Newtonian potential in which case the vDVZ discontinuity does not appear. Therefore, vDVZ 
discontinuity does not appear in a curved background, since we can compare the smallness of the 
mass of graviton with another scale that is the cosmological constant. We can also generalize 
the Newtonian potentials to D dimensions. For this we need the general Green's function for the 



operators ^ and 



92- 



They are 



G (x, x) 



(2n 



(D-5) 

2^^r 



D 



and 



G (x, x') 



1 



1 



D-l D-3 

(27r) 2 r 2 



1 



3-D 
4 



K d-3 {rrnf^ 



amplitude equation ( 85 ) as 



respectively. T is the gamma function and Ky is the modified Bessel function. We first write the 

^5{x - X2) 



(Z?-2^Kmim2 . , 

d[X XI j ^2 
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mt, 



{D-l){D-2) - ^^^5^^^^ - 



K mi m2 
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1 



{D-l){D-2) 
Taking again static sources, we find the potential as 



Q2 _ 



6{x - X2). 



(89) 
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(90) 



For any desired dimension this equation gives the Newtonian potential in flat spacetime and without 
any mass term. For k > the first term signals the ghost problem for all dimensions and the last 
two terms give attractive forces. We can confirm j21| that in any dimension the quadratic curvature 
gravity theory is not unitary except for the NMG point. 

For curved backgrounds the calculation of the Newtonian potential is complicated because of 
the Lichnerowicz operator. Nonetheless, we can say that in curved backgrounds there will be an 
additional term to the Newtonian potentials. Apart from the term that contains Lichnerowicz 
operator, other terms can be calculated as above. 



XI. CONCLUSION AND DISCUSSION FOR CHAPTER-2 



The tree level scattering amplitude of the general quadratic gravity theory with a Pauli-Fierz 
mass term in D dimensions is found and various limits of this amplitude equation is studied. For 
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these limits the Newtonian potentials for three and four dimensions are calculated. First the flat 
space limit is taken and then the massless limit. Looking at the potential energies it is seen that 
the theory is not unitary. For taking care of the ghost term the /3 — t- limit is taken. In this limit 
the theory becomes unitary but the theory suffers from the vDVZ discontinuity in four dimensions. 
Massive gravity gives the correct Newtonian potentials in three dimensions since the M"^ ^0 Hmit 
does not exist. The second limit that is taken is the massless limit and then the flat space limit. 
Analyzing the pole structure and the residues of the amplitude and calculating the Newtonian 
potential energies it is seen that the theory is unitary for a special condition for the coupling 
constants of the theory that is 8a + 3/3 = in three dimensions. Also, the coupling constant 
of Hilbert-Einstein term becomes negative. Apart from three dimensions the general quadratic 
theory is non-unitary. To make the theory unitary again the /3 — )■ limit is taken and for this case 
it becomes unitary. The usual Newtonian potential is obtained in four dimensions by taking all 
coupling constants of higher derivative terms to zero. Also, the Newtonian potentials for generic 
dimensions in flat space time are calculated. 

Part III 

CANONICAL STRUCTURE OF HIGHER 
DERIVATIVE GRAVITY IN 3D 

XII. INTRODUCTION 

In the previous chapter we saw that the general quadratic gravity theory has a special case in 
three dimensions. This special case is 8a+3f3 = and k < for the Lagrangian k 
and named as New Massive Gravity (NMG) |35| Also, the parity violating extension of NMG 
is found j35| 155] by adding a gravitational Chern-Simons term to the theory. These theories have 
massive ghost-free spin-2 particles in their free spectrum around both flat and (anti)-de Sitter 
(A)dS. The most interesting property of NMG is that it gives a non linear extension to the Pauli- 
Fierz mass term for spin-2 particles. It is the only known example which provides this nonlinear 
extension. Thus, NMG is a suitable candidate to be a perturbatively well defined quantum gravity 
theory in three dimensions if it is unitary beyond tree level. 

The theory is studied in many directions. Its ghost-freedom and tree level unitarity |17| [27| l35l 
1351 E2] and Newtonian limits [27 have been studied out. Furthermore, its classical solutions and 
issues related to the classical solutions were studied p5ti38| H0tii2] . Also, its supergravity extension 
were given in |75j. 

In this chapter we study the canonical structure of the general quadratic curvature theory and 
give an explicitly gauge invariant analysis of it in three dimensions. The analysis is done for both 
flat and de Sitter (dS) spacetimes. In the flat space analysis the gravitational Chern-Simons term 
is added to the theory. By writing the canonical form of the general quadratic action we can easily 
see how NMG is singled out among all other higher curvature theories as a regular "harmonic 
oscillator" which can be thought as massive free field. These type oscillators do not have the 
Ostragradskian instability which spoil every higher time derivative theory [76J^'^. Except of NMG, 
at the linearized level, all the other quadratic theories are ghost-ridden Pais-Uhlenbeck oscillators. 
Moreover, the discussion is extended by discussing the Newtonian limits, weak fields, the scattering 
of particles with mass and spin at the tree level. 



It was claimed that adding interactions might yield stable higher-time derivative theories |77| . 
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The content of this chapter is as follows: In the first section the flat spacetime analysis is 
done by writing the canonical structure for both quadratic curvature theory with and without 
the gravitational Chern-Simons term. Also, the effects of static sources are analyzed in detail. 
Moreover, the weak field solutions are obtained in circularly symmetric case. In the second section 
the discussion is extended to curved spacetimes namely to de Sitter spacetime. Some of the useful 
calculations are given in the Appendix. 



XIII. HIGHER-DERIVATIVE SPIN-2 IN FLAT SPACETIME 



We start the analysis with the flat space case since it is simpler than the curved space background 
which is discussed in the next section. Without any constraints on the parameters the general action 
is 



-R + aR^ + I3R^ 



(91) 



To get the desired spin-2 model the action must be expanded around flat background gfj_i, = 
'HpLv + h^y^ where r]^^ is the usual flat spacetime metric with signature (—,+,+) and h^y is the 
spin-2 field which actually is a symmetric rank-2 tensor and without any constraints has not only 
spin-2 component but also spin-1 and spin-0 components. In this section we do not consider the 



gravitational Chern-Simons term. However, in the next section we add it to (91). To get the 
action for h^y we first linearize the field equations that comes from the variation of (91 ) and then 



integrate the linear field equations. While integrating them, one must pay attention to the overall 
sign factor since the signs will become important in the discussion of unitarity. Up to boundary 



terms (91) becomes 



^ I d X h^y 



df^d") Rl + pagi^ 



(92) 



Here, Q'^'^ is the linearized Einstein tensor and Rl is the linearized Ricci scalar which is defined as 
Rl = {gfj,yR^'^)^. The linearized Einstein and Ricci tensors, and linearized curvature scalar are as 
follows 



R 



'l' - ^V^^Rl, n.L = UaUp 



Rl = dadgh''^ - ah, 



h = 7?^" V> 



(93) 



where □ is the D'Alembertian operator that is □ = d^d'^ = —d^+V'^. Since the metric is perturbed 
around fiat background, we raise and lower the indices with the fiat metric Tj^y. To analyze the 
canonical structure and explore the free fields h^y must be decomposed. This is done as follows: 



/loi = -(^ijdjT] + diNL, 



h. 



00 



(94) 



where (j), x> ^^Vi ^L, N are free functions of time and space {t, x) and di = V^. 

After decomposing the spin-2 field, the components of the linearized Einstein tensor must be 
written in terms of these free functions. Let us give some details of these calculations. 

We first write the trace of the spin-2 field: 



(95) 
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Putting ( 94 ) in ( 95 ) gives us 



-N + h\ = -N + 4) + x, 



where we have used did^ 



-1 and note that K. 



can be written as 



- X- The hnear Einstein tensor has three com- 



ponents that are Qqq, Qq^ and Qj^:. Using (93) and summing the repeated indices these components 



After using ( 94 1 these components can be written in terms of gauge invariant objects as follows 

1. 



^00 



^ ^(Sij + didj^ q - didjc, 



where " ' " means differentiation with respect to time t and g, o", and (p are defined as 

q = V^N - 2V^Nl + X, cj = i- V^rj, 



(96) 



(97) 



which are gauge invariant under the transformation (5^ /i^i/ = d^C,u+duC,^ and the Laplacian operator 
is = df = did'^ = didi . Note that, cf) is already a gauge invariant component of h^y. 

Let us show this invariance of these objects: First we define the three vector as = 
(Aq, Ujdj + diK) and then we transform the components of the spin-2 fields. From the trans- 
formations of /loo c^nd /loi we get 



7?' = r/-A, N'L = NL + k + Ko, N' = N + 2Ao 



(98) 



From the last component we get three more equations by multiplying h[j = hij + di^j + djS^i with 
djdi, Eiidi and 5ij that are 



(99) 



respectively. From here we can see explicitly the gauge invariance of (j). Using ( 98 ) and ( 99 ) , one 



can get (97). The linearized Ricci scalar can be written in terms of the gauge invariant functions: 

= dohi + df (/loo - V) - '^dodihoi + didjhij, 



again using (94) the scalar curvature becomes 

7L 



+ X - V2(/) + V^iV - 2V^Nl. 
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Putting the first equation of (97) finally we get 



(100) 



Although we have six free functions, the Bianchi identity d^Q^ = reduces the number of arbitrary 
functions to three. Therefore, there are two sets of functions that can be used which are u, q 
and (/>, cr, Rl. We can write the total action in terms of these gauge invariant free functions. Let 
us start with the Einstein-Hilbert part of the action that is 



Ieh 



1 

'2k 



then doing the summation about the repeated indices and using (94), (96) and (97) we get 



Ieh 



(101) 



Here we can see that the pure Einstein theory does not have any propagating degrees of freedom 
since there is not a dynamical part which shows itself as a D'Alembertian operator. With the same 
route we can find the other two parts of the action. Instead of using this technique one can also 
use the self-adjointness of the operators. In both ways one can write the quadratic parts of the 
action in explicitly gauge invariant forms. They are 



2a + 13 



2a + (3 



d^xRl, 



(102) 



d^x h 



d^x 



+ -Rl 

2 ^ 



J d^x {qa(P + aaa) . 



(103) 



In the first action we move all the derivatives on h^^, and that combination gives us again the 
linearized Ricci scalar. In the second action /«, first the derivatives are moved on the spin-2 field 



and then (93) is used. Finally using the Bianchi identity one can achieve (103). Adding (101), 



(102) and (103) the total action comes out in gauge invariant combinations as 



d-^x 



+ (2a + /3) {q - n(f,y + Pqa^ 



+ 



d-^x ( aUa H -a' 

up 



(104) 



From this equation, it can be seen immediately that a shows itself as a single massive scalar field 
--'o. Not to have negative mass k/3 must be negative and not to have ghost (3 



with mass 



must be positive for the a field. For these signs k must be chosen negative. For the remaining 
part of the action the discussion bifurcates whether 2a + ;8 = 0, or not. Let us discuss these cases 
separately: 



A. 2a + (3 case: 



In this case the q field can be eliminated by taking the variation with respect to it. Doing so 
gives 



1 



2{2a + p) \k 



+ /3n0 +□ 



(105) 
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and inserting the q field into (104), the action for the (j) field becomes 



(fx 



/3(8a + 3/3) 



+ 



(4a + 13) 



(106) 



4(2a + /3) ' ' 2K{2a + py ~^ 4^2 (2a + ^)^ 

In this action there are some special limits. The first limit is the NMG limit 8a + 3/3 = 0. For 
this limit the higher derivative term drops out and the (j) field describes a single massive degree of 
freedom 



NMG,<f> 



(107) 



A/a («+;i;. 

which has the same mass as the a field. Since both degrees of freedom have the same masses in 
the NMG limit, the theory is parity-invariant. Note that, not to have a ghost k must again be 
negative. Also, NMG is not a higher derivative theory at the linearized level. The theory does not 
have the Ostragradski instability which appears in theories that have higher order derivatives. 
Another interesting limit is 4a + /3 = 0. The middle term in (106) drops out and the action 



(106) becomes 



This theory is tachyonic for «; < and gives a ghost term for k > 0, keeping in mind that k/3 < 0. 
Taking /3 = also drops the higher derivative term. For this case the action reads 

U = — (fx ( ^—Ci 



(108) 



(109) 



In this case, to avoid ghost term k > and not to have tachyon a > 0, where the mass can be 
defined as = 



For the general coupling constants, the above action (106) is a higher derivative Pais-Uhlenbeck 



oscillator. In order to decouple the fields we define new fields that are linear in terms of old fields. 
By inspection the fields can be defined as follows 



□c 



^2 



□c 



ni't 



(110) 



and with these definitions the uncoupled Lagrangian must be as follows 



Kiipi (u - mfj + K2ip2 (□ 



Here Ki and K2 are unknowns. To find these unknown terms we put (110) into above Lagrangian 



and compare that result with ( 106[). A fter getting the unknown factors the action can be written 
in terms of simple oscillators and (106) becomes 



1 

64k {2a + /3y 



d^x 



(8a + 3/5)2 (^^ _ ^2^ _ p2^^ (^^ _ ^2^ ^^j ^ (^^^^ 



and the two fields (pi, (p2 have masses m<j and rUg respectively. rUg is defined above and is 
defined as 



1 



K (8a + 3/3) 

Note that for the general case the only restriction on the coupling constants is 2a + /3 7^ 0. For 
K < the (fi field gives ghost and for k > the (/?2 has negative kinetic energy. Also, not to have 
negative mass term 8a + 3/3 must be positive for positive k and negative for negative k. 
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B. 2a + /3 = case: 



In this case we can follow two routes: we either start from (104) or from (111). These two 



actions give different results and as we have seen above in the second action this limit provides a 
singular theory, w hile i n the first action it does not. 



We start with (111) and take e = 2a + limit. By this limit 



where in the last step we have taken the Taylor expansion around small e. Note that Tni does 



not depend on e. With the same limit we can write the decoupled fields in terms of e as ip2 
(b „ P*^ . X PS + ir^ and again does not depend on e. With these expansions, (111) gives 

us up to second order 



1 

8Ke 



(n - mfj (/) ^ - 4e0 (□ - mfj (p + O (e^) 



(112) 



This action is a degenerate Pais-Uhlenbeck oscillator which is known as ghost free for some param- 

(113) 



eter ranges. However, if we start from (|104|) and set 2a + /3 = 0, we get 



^ ' d^x {qUcf) - TTi^gC^ 



Taking variation with respect to g or gives us a massive wave equation. However, these equations 
do not say anything about the ghost structure or tachyonic behaviour of the theory. So that, we 



must again separate the fields by redefining the fields such that q 



m' 



(^1 + ^2) and = -^2- 



With these definitions (113) turns into 



d^x 



(^1 



(^5 



(114) 



From the above discussion we know that (3 must be positive not to have ghost term, but in this case 
if it is taken positive, ^2 becomes a ghost excitation. As we discussed in the previous chapter the 
Newtonian limit of this theory is interesting [27] . The Newtonian potential of this theory becomes 
zero when two static sources are taken into account, since the ghost excitation that is the repulsive 
component cancels the spin-2 part which is the attractive component. This is the same situation 
that happens in the pure Einstein gravity. 

C. Adding static and spinning sources 

Up to now the analysis depended on theories without any interaction. However, when interaction 
enters into the picture it may change the particle spectrum of the theory. So that, we turn our 
analysis to source dependent higher derivative gravity. First we add static sources to our analysis 
and then we generalize this analysis by adding spinning masses in flat spacetime. 

1. Static Sources: 

The matter can be added to the theory by the usual gravity-matter coupling that is 



- I d'^xh^.T'"'. 



(115) 
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For the static source case we take = p (x), = 0, T*^ = 0, and the above action becomes 



(116) 



1 /■ o 1 

/source = - / d xNp{x) = - 

where we used q = V^N — + x ^ N = + 2Nl — ^X- After dropping the boundary 

terms and taking the source terms static, that is p (x) = 0, and using the symmetry property of 



the Green's function in the first term of (116), the source action becomes 

-^source — ^ 



d xq^p. 



We add this action to the general total action (104) 
1 
2 



d-^x 



+ {2a + /3) (g - □0)^ + /3gn0 + q^p 



(117) 



where is the action that is only constructed from the a field and is already in the simple harmonic 



oscillator form. To eliminate the last term in the parenthesis of (117) we define new fields so that 



we can decouple the fields easily. The last and the first terms can be combined and redefined as 
a new field (p = (f) + K^p and by inspection the q field can be redefined as q = q + up. Putting 
these new fields into (117) and taking □ = —df + V^, we get the total action as 

^ -(ipq-Kipp + a'^)+{2a + l3){q-Uipf 



d^x 



-/? {^qOip — npOip 



Kqp + K^p^ + aOa 



Let us concentrate on the NMG case, that is 8a + 3/3 = 0. For this case 2a + P 
variation with respect to q field, the following equation comes out 



(118) 
Taking 



-if + ^ {q -n^) + (3 {Oip - Kp) 



0. 



From this equation q can be defined in terms of the other fields as 

2 



q = 2kp — Uip 



K/3 



ip. 



Putting (119) into (118) the full action becomes 

1 
2 



d'^x 



:ii9) 



(120) 



The last term in the action is the interaction part and it gives an attractive potential energy for 



negative k. In order to find the potential we take the variation of ( 120 ) with respect to (p field that 
gives 



Putting the ip into the interaction part of (120) the potential energy reads 



(121) 



U 



d^x pi 



V2 



8tt 



(122) 



From this equation it can be seen easily that for negative k the interaction part gives an attrac- 
tive potential energy, that means in NM G cas e adding static sources does not spoil the unitarity 
structure or the tachyonic behaviour. In (122) the sources are defined as pi (x) = miS^'^^ {x — xi), 
P2 (x) = m2<5*-^^ (x — X2), and Kq is the modified Bessel function of the second kind. This result 
also matches with the result that we found in the previous chapter [27] . 
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2. Spinning masses: 



For this case the energy-momentum tensor must be written as 



Too = ■m.^*-^'* {r 



n), 



0, 



where m is the mass and j is the spin of the point-hke source. The general amphtude and Newtonian 
potential was calculated for static sources at flat spacetime in Z)-dimensions in the previous chapter 
||27j . They become 



AA 



-2T' 



rpflU ^ rj-i/ 



K 



T 



-T' 



1 -1 



ia + 3/3) 9^ - -d^ 

K 



and 



U = ^ ^ [Ko [r mg) - Kq (r ma)\ 



in three dimensions. The only added part from the spin will come from the Tqj components of the 
energy-momentum tensor. It will read as 



41^0 ( + -5' 



niO 



pm'^ 



r — ri I 



1 



Q2 Q2 
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r - r2 1 . 



"9 \^ ^ "9/ 

The Newtonian potential of these operators was calculated in the previous chapter. Using that 
result the potential energy reads 

4^/ 1 



iO 
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Jlj2 



52 



Using e^^e^^ = 6^^ the spinning part of the amplitude becomes 

T 



Pm'g 



,2 1 ^iO _ JlJ|5^52(^-._^-.^)_^^^^2(^-._^-.^) 
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J1J2 a j:2/- - n 



(9^ — 



di6^{r- r2). 



Using the Green's function of these operators and taking the differentiations in front of the functions 
and carrying the space integrations, we get 



J1J2 ^2, I- - I 
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Note that V'^ln\fi — = — 27r(5^ (n — and if we assume that the sources are in separate 
positions, that is fi ^ r2 the first term in the above equation vanishes and the second term becomes 
{y'^ -mfjKQ{mg\ri-r2\) = 27r(5^ (n - = V^Kq ("igln - ^2!) = m^Ko ("igin - ^21) 
and the amphtude takes the form 

- 47^0 (^a^ + ^52) ' = -||Ko (m, |fi - f2|) . (123) 

With this spinning part, the total Newtonian potential energy, [/ = Amplitude/time, can be written 

as 



^ = ^ (mim2 + 4mgjij2) Kq {nig |fi - r2|) - ■^mim2Ko {rus \ fi - f2\) . (124) 

The signs of ji and j2 are not fixed so that they can be both positive or negative. Therefore, 
potential energy coming from this part could be both attractive or repulsive. In the NMG case 



this situation does not change since for this case only the last term in (124) drops out. 



D. Weak field approximation 



Up to this point we have found some results from the linearized theory. In this section we will 



try to get some of the above results from the nonlinear theory (91 1 . We start with the ansatz 



ds^ = -f (r) dt^ + Kv^dr"^ + r^dO^, 
f{r) 



(125) 



and insert it to (91). Here / (r) and 6(r) are functions and the action will be varied with respect 



to these functions . We will consider only the NMG theory. To find an approximate solution the 
functions can be defined as follows: / (r) = 1 + dra{r), 6(r) = 1 + drv{r). Here, v (r) and 
a{r) are small functions. Putting the ansatz into (91) with the defined functions we get 



-V + 2pv" + 2/3a" + rpa'" = 0, 



(126) 



fir'^a" + -r^a + 2rl3v' - 2(3v = 0, (127) 

K 



up to first order in v {r) and a(r). Here, ' denotes differentiation with respect to r. From (126) 



and (127), v (r) can be written in terms of a(r) as v (r) = a{r) + £a'(r). To have an ordinary 



differential equation we put v (r) into (127) and we get 



r^a" + ra - a (mlr'^ + l) = 0. (128) 

By solving this equation we can find a (r) as a (r) = ci/i {rrigr) + C2K1 (mgr), where ci and C2 are 
constants. From this equation also v (r) can be determined using v (r) = a (r) + |a' (r). From the 
above ansatz we know that (700 ~ ~1 ~ dra{r). For we first take the square of b{r) and 
drop the second order term. Then we expand up to first order that is ~ 1 — e. Combining 
these approximations we get grr ~ 1 + dr [2v (r) — a (r)] . We know that for decaying fields that 



See the details of this Weyl trick in [7S] 
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is for r — >■ oo, a (r) — )■ 0. Since Ii diverges to infinity for decaying fields ci must vanish, and the 
metric components become 



500 



-1 + cKo {rugr) , ^1+ dKi {nigr) 



(129) 



Here c and d are constants. When we compare these components with Schwarzschild solution we 
see that the constants must be related to the mass of the sources. This is consistent with our 



earlier result (122). 



E. Higher-derivative gravity plus a Chern-Simons term 

In this part the above discussion will be extended by adding the gravitational Chern-Simons 
term to the quadratic gravity theory in flat space. We will follow the same route; first we will 
decompose the linear form of the action that comes from the gravitational Chern-Simons term. 
Then the coupled field will be decomposed by taking the Fourier transform of the Lagrangian 
and writing it in a matrix form and diagonalizing this matrix. The general quadratic action with 
gravitational Chern-Simons term is 



/ 



-g 



R + aR^ + PR 



1/2 



2/i 



where eoi2 = and /x is the Chern-Simons coupling with an arbitrary sign^^. 
The linearization of the gravitational Chern-Simons bit yields 



Ics 



1 

'2^ 



d'xe^a.^G2''d^h''^ + 0[h^). 



First we write all the summations that is 



(130) 



= ^ y c^^a; dj {^-Qiodohjo + Goodikjo - Qj^dihoo - Qj^dfihik + Gokdjhik - G^k^jhok) , (131) 



where eij = eoij . Using ( 94 ) and ( 96 ) the terms that appear in the last line in the above equation 
are 



— eijGojdihoo — - 



^ijG^kdjhrk = 2 [(yV'(t) + iV' 



(132) 



where we have taken out total derivatives in the needed steps and dropped the boundary terms. 
With these equations the action of the gravitational Chern-Simons term takes the form 



^CS = ^J d^x [a {q + 



(133) 



Without the a, /3 terms, but with a Pauh-Fierz mass term, canonical analysis was carried out in [591 168| . 
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where we have defined the q and a terms in the previous section. With these, action ( 130 ) becomes 
. 1 



dfx 



1 



((t)q + CT^) + (2a + I3){q- Ucjif + p {qUcj) + aUa) + -a{q+n 



(134) 



where all the fields are gauge invariant. For the general case, that is 2a + /3 ^ and both 
a 7^ 0, /3 7^ 0, we eliminate the q field by taking the variation of (134) with respect to the q field. 
From this equation of motion we get 



and 



= -(A + 2 (2a + /3) (g - net)) + pUff) + ^a, 



ia + P)U(j)- -(f>- -a. 



(135) 



(136) 



Putting (136) in (134), the general action takes the following form 



(fx \p 



1 1 



1 1 

/3k (8a + 3/3) 



+ 



1^ (4a + /3) 



4(2a + /3) 



// 2/x(2a + /3) 
(4a + P) 



2{2a + py 



4K(2a + /3) 



(137) 



For a proper analysis we must decouple the a and the (j) fields. This decoupling procedure can be 
done mathematically for generic a and (3 but since we are interested in the NMG case we take the 
8a + 3/3 = point. For this case the action reads 



Inmg- 



cs 



(Tx 



1 



^2/32 



a 



+ 



2?TZq 9 / 9 ' 

a(j) + nig {<j)0(p — nig^) 



To decouple these fields we will follow a different route: First we take the Fourier transform of the 
fields and put the Lagrangian in a matrix form then we diagonalize this matrix. 

Taking the Fourier transform of the Lagrangian and putting in the matrix form yields 



Lft = ( 5" rugcf) ) 



lltg 




Using the eigenvalue equation det {XI — A) = the eigenvalues of the A matrix can be found where 
/ is 2 X 2 identity matrix, A is the above matrix and A is the eigenvalues of the A matrix. Solving 
the eigenvalue equation for A we get 



+ 



1 



/i2/32 



X+(k 



2 , 2 



+ 



0, 



with roots 



A± = -k^ - 

Now, let's find eigenvectors. For A-t; 



2/i2/32 ^ 



I rrig + 



2«2- 



(138) 
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N± 



ml + 



, N± 



1 + 



mn 



(140) 

where N± are the normahzation factors. Then we construct a modal matrix from the components 
of the eigenvectors such that 



P 



( ' 



2 Jim, mg \J'^l + ij^) 



(141) 



and since the eigenvectors are orthogonal to each other the inverse of the modal matrix is equal to 
its transpose. With the inverse of this modal matrix we can define new fields such that = P^^X, 
where X is the original field, 



iV_ iV_ 



2nl3mg mg\J"''9 ' 4^t^ 



mg4> 



(142) 



Then, the diagonalized A matrix can be constructed as D = P ^ AP which is just D 



A+ 
A_ 

from which the transformed Lagrangian reads LpT = ^ D^. Then, the action in the decoupled 
form becomes 



Inmg- 



■cs 



or in a more concrete form 
Inmg-cs 
where the masses read 



m„ + 



1 



1 



1 



1 



(143) 



(144) 



(145) 



and the fields are inverse Fourier transformed. The masses are same with those of j35| [551 E5]- 
Since masses of the helicity modes are different this theory is a parity violating theory. To check 
the results of topologically massive gravity we take ^ — )■ limit. In this limit m+ diverges and 
drops out, therefore we are left with a single degree of freedom that has a mass m_ = — / n 
3j. This result can be seen by Taylor expanding the square root part of (145) up to third 



order. 



XIV. HIGHER-DERIVATIVE SPIN-2 IN A DE SITTER BACKGROUND 

Up to now our analysis was based on the flat space time. Now, we will change the background 
to a constant curvature background. Specifically we will study the canonical structure of higher 
derivative gravity whose action is defined as 
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1 



{R - 2Ao) + aR^ + 



(146) 



in an (anti)-de Sitter background. Here Aq is the bare cosmological constant. The Hnearization of 



( 146D yields 

1 
'2 



-gh 



+ (2a + (3) ^^^D - V^V^ + -^g'^'' ] Rl + f3 { 



-,J1U 



£2 



(147) 

where 1/^^ jg the cosmological constant and a = ^ + j^a + ^(3. The cosmological constant 
can be related to the coupling constants a, (3, k and the bare cosmological constant as = 

^^^3^^^^ 1 ± y/1 - 8kAo (3a + /3) [271 [721 |T2]. We will carry the analysis in a de Sitter (dS) 
background since it is simpler than an anti-de Sitter (AdS) background. Nevertheless, from these 
expressions one can find the results in AdS spacetime by taking i ^ ii transformation. This 
transformation can be taken since the results that come in the dS background are analytic in £^^. 
For dS, we define the background metric cjfj^u in the Poincare form as 



ds' 



e_ 

t2 



-dt^ + dx^ + dy^ 



(148) 



by which all the raising and lowering operations and covariant derivatives will be made. The 
perturbation can be defined as 



(149) 



where jiVfj^i^ — dtJ-f ^^^cl Tj^y is the flat spacetime metric. With this perturbation the linearized 
forms of Einstein and Ricci tensors and scalar curvature can be written as 



R 



'111/ "flU 



Rl - ^ V' 



R 



(yV^.Ka + V"V^ V - ° V - V;.V^/i) , Rl = VaV^h""^ -ah- j^h, (150) 



where D'Alembertian operator defined as □ = V^V^ = j^ri'^'^'V f^V^. The deviation part of the 
metric h^i, can be decomposed into its "spatial" tensor hij, vector hoi and "scalar" /iqo parts as 
follows 



t2 

fj2 



riV, 



(151) 



To keep the signature intact, one also needs to Wick rotate a space coordinate 
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where Vj = V^/ y^vf and the covariant derivative is defined for the two-dimensional space metric 

7ij = ^Sij. The Latin indices are i = 1,2 for space dimensions. Since the components of the two 

dimensional metric has no space dependence and is flat, the covariant derivative reduces to partial 

derivative, Vj — )■ di, and di = di/'^—d'^. iik is the Levi-Civita tensor which is related to the 
corresponding tensor density as 



t2 



(152) 



The convention for e^fc is ei2 = 1 (the convention for Levi-Civita tensor density for the upper indices 
is e^^ = 1 naturally with the induced metric). Therefore, the final result of the above decomposition 
becomes 



t2 



(dij + didj^ 4> - didjx + (eikdkdj + ejkdkdij ^ 



hoi = ^ i-eijdjr) + diNi,) , ^oo = :^N. (153) 

The raising and lowering operations for spatial indices are done with 6ij. Note that, the specific 
choice of decomposition involves |^ coefficients and these coefficients help us to check the flat space 
limit, £ — >■ oo, I ^ 1 at every step in our calculations. 

At this stage we will find the gauge invariant combinations that will be constructed from the 
six scalar functions. In fiat space case is gauge invariant. However, for curved background it is 
not gauge invariant anymore. The components of hfj,i, transforms under the gauge transformations 
Sch^i. = Vi^Cu + ^uCn as 



2^ 



f2 



£2 



K + Co + 



5cN 



Co + -Co 



(154) 



where the components of are defined as C/i = {C,Oi~eijdjC, + din). Looking at the linearized 
Bianchi identity, V ^Q^^ = 0, there should be at least three independent gauge invariant combina- 
tions which are constructed out of the six scalar fields and their derivatives. These combinations 
can be found by inspections. However, looking at the independent components of the gauge- 
invariant tensor the gauge-invariant combinations can be found easily. Following this route 
four gauge-invariant functions can be found: 



2.. 11 



p=-[(t>--N 



N ■ 



a = 



and from the Bianchi identity we get a relation between these functions as 



p-q = 0. 



(155) 



(156) 



We can find the components of the linearized Einstein tensor in terms of these gauge invariant 
fields, that are 



='00 ■ 
t 

'21 



eikdkcr) 



2£ 21 
(Sij + didj^ q - didjp - (eikdkdj + ejkdkdi^ a 



(157) 
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We can also write the linearized curvature scalar in terms of the gauge-invariant fields 



Rl 



£3 



(f-p) 



(158) 



In the second equality we used the Bianchi identity (156). 



Using the above relations (147) can be written in terms of the gauge-invariant functions. The 



Einstein-Hilbert part of the action takes the following form: 



Ieh 



d X \J-g h^yQ^ 



fit/ 



t 



t2 



fRL + -[fV'f + p'+a^ 



(159) 



For the 2a + (3 part of the action, using the self-adjointness of the involved operators, h^i, can be 
replaced by some gauge-invariant combinations. Also, doing so makes the computations simpler 
as in the flat space case. With this trick, the 2a + f3 part of the action becomes 



(2a + /3) 



qKu [ g^'a - V^V^ + ] Rl 



(2a + P) 



d'^x 



-gRl 



L- 

(160) 



For the /3 part of the general action, we have 

1 



d-'x 



-,au 



r^RL 



d^X 



;hRL 



First, writing i?^^ (150) such that the indices of covariant derivatives /i and v stay at the left, and 
then using the Bianchi identity, the following action can be found for the (3 part 



d-^x 



-g 



-K^g'l'' + \rI + Ik^q'l 



If we had not used this method and computed h^yUQ'^ directly, to put the result into an explicitly 
gauge-invariant form would have been more difficult and taken more time. Finally, the general 
action can be written in terms of the gauge-invariant form by collecting all the parts that are 
computed above as follows 



1 



d-^x 



a + 



2_l 



f t 



+ aVV + p2 + pV'^p + (VV) 
RlP - pV^f 



+ (2a + P) ^Rl 



(161) 



f3 f3 



The above action gives ( 104 ) when the flat space limit is taken. Also, the fields that appear in ( 161 ) 
are not independent. However, defining new field such that = V^/ , and after using the Bianchi 
identity (156), the above action (161 ) takes a rather simple form in which the a field decouples 

f5 



1 



d-'x 



2/3\ t 



{p^ - pV'^p - Lp^ - ttpV'^p - tptp - (/^p) I -h la 



(162) 
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where the a action is 



j d^x (d^ + (tVV) + {a + 



21 



(163) 



In flat spacetime, cosmological Einstein theory does not have any propagating degrees of freedom 
for vanishing a and /3 in D = 3. Just hke the flat spacetime, without higher curvature terms the 



total action does not have any propagating degrees of freedom. For non- vanishing a and /3, (162) 



has three degrees of freedom. From the decoupled fleld we can read the mass if (163) is put in a 



correct canonical form. For a minimally coupled scalar fleld the correct canonical form is 





f d'x < 











+ [di'^y 



(164) 



To get the correct canonical dimension the a fleld is rescaled as a 



£2 



a. After this rescaling, 



(163) becomes 



t3 



a 2 

;5 ^ F 



a 



and for the a fleld the mass reads 



2 

£2 



1 



12a 



4 



1 



4(3q-/3) 

w ■ 



(165) 



(166) 



Unlike the flat space case, diagonalizing the ip, p action is complicated for the general coupling 
constants. Since we want to see the oscillators of this theory we can use other methods. One of the 
methods is to Fourier transform the flelds in the x space and then computing the zero-momentum 
limit, that is dropping the Laplacians in the action. The action does not change the number 
of degrees of freedom by doing this manipulation, since, (fleld) is not the lowest order term 
anymore. Another method is to directly work on the equations of motion. Both methods will be 
studied separately below. 



Masses from the nonrelativistic limit 



In this part we take the nonrelativistic limit of the general action ( 162 ) by dropping the V 



terms. Since the a part of the action is already in an oscillator form, it is not taken into account. 
Then the action becomes 



/ = ^ I d^x 



(167) 



At this step we must decouple the (p and p flelds. For this reason we rescale the fleld as ip ^ip. 



Also, we separate the constant 2a + /S = f + ^^±^ to see clearly the NMG Hmit. With these (167) 



becomes 



d^x 



4P 



•2 _ ^ 
t2 



(8a + 3/3) / .2^ 3v92 
+ 4 +1^ 



(168) 



42 



However, this rescaling is not enough to separate the fields. Therefore, we define a new field 
$ = — 2p. With this definition the if and ^ fields decouple. The $ field action can be read as 



(Fx 



$2 _^ 



t 



a 2 



(169) 



and gives the same mass as the a field (166). Therefore, the $ field becomes the spin-2 helicity 



partner of the a field. The action for the field reads 



1 

ia + 3/3) 



t ( 24a 

a ^ 



6/3 \ 2 



(170) 



which is the spin-0 mode. To read the mass of this mode, we must put it into the canonical form. 
For this reason we again rescale the (p field — >■ After this rescaling the mass can be read as 



1 



4 / 3a + /3 



AC (8a + 3/3) 



8a + 3/3 



(171) 



For the NMG case this mode drops out and we are left with the rrig and this result matches with 
[35j. Note that in [35^ the analysis was carried out by introducing auxiliary fields, and these fields 
can be eliminated in such a way that the action gives spin-2 field with a Pauli-Fierz mass. However, 
we reach the same results by using canonical analysis. For general coupling constants, to make the 



same analysis as in [35j we must introduce two auxiliary fields and rewrite the action (91 ) in terms 
of these fields: 



1 



R- rGau-(t^R + 



+ 



(172) 



where (j) and f^^^ are auxiliary fields and the masses are 



linearize (172) around flat background, that is 



' k(8cH-3/3) "^2 



- \ . Then we 



linearized 



4>Rl 



K (8a + 3/3) 



(173) 



Here we can see explicitly that for the NMG limit the (j) field drops out and /^jy can be chosen as 



-h^u- However, for the generic case it is not clear how the fields (/>, f^u and h^y decouples. 
One possible way is to rescale the h^u field [80j . 

In this section we have discussed the canonical structure of the generic action in the nonrela- 
tivistic limit. It would also be interesting to get the same results from the analysis of the relativistic 
equations. The next section is devoted to the discussion of the relativistic equations for the NMG 
limit. 



B. Equations of motions in the NMG case 



The action (162) for the 8a + 3/3 = case becomes 
/3 /■ { ot 



I 



+ 
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where we dropped the a field. From this action one can claim that there are two degrees of 
freedom, which is in conflict with our earlier results \ib\ 155] . However, if we look at the Hessian 
matrix?^ = gg-, that is 



4P 



-2t 



-2t 4 



(175) 



and compute its determinant, defH = 0, we see that there should be a constraint in this model. 
Therefore the time derivatives of the fields cannot be separated in terms of the canonical momenta 



a^ = 4Fl^-^^-t^, 



Since the equations of motion are needed we take the variations of (174) with respect to and p 
fields which yield 



5ip : 



V 



{2lp + tV' 



p + p 
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2P 
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6p : 



{t^ - 2p) 



2P 



[if - V^p) - 2t^p 



rdo {2p - tip - ^p) 



0, 



0. 



(176) 



:i77) 



To decouple the fields, we define ip = V'^p with a hint that this choice makes Rl = 0. With this 
definition, the other equation becomes 



-ip- 



V9 + vv 



t3 



m„ 



1 



ip = 0. 



(178) 



This equation is not in the canonical wave equation form in dS. We again rescale the field p ^ ip/t 
to have canonical wave equation and this rescaling yields 



ip+^ip + VV) - ^^rn^ip = 0,^ (a-mfjp = 0, 



(179) 



which is the same as the a field. From the relativistic equations, we again see that for NMG limit 
the quadratic gravity theory is parity preserving theory in three dimensions. 



XV. CONCLUSIONS AND DISCUSSION FOR CHAPTER-3 

In this chapter the canonical structure of the linearized quadratic gravity theories have been 
studied. The analysis is done in an explicitly gauge-invariant way. Moreover, the analysis is 
made both for flat and dS backgrounds in three dimensions. In fiat background case, the general 
quadratic action has been written in canonical wave equation form. After the fields are decoupled 
they generate three harmonic oscillators. The coupling constants k, a, (3 are chosen to have a 
ghost-free and non-tachyonic theory. When the coupling constants are fixed in this way, the NMG 
theory is singled out as the unique unitary higher derivative massive gravity (but not a higher- 
time derivative theory). Apart from this theory, all other higher derivative gravity models are 
higher- derivative Pais-Uhlenbeck oscillators which have ghost modes. 

The analysis is also extended by adding static sources and spinning masses to the theory. 
The effects of the sources are described by computing the Newtonian potentials for both cases. 
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In addition, the weak field limit of the theory is computed at the nonlinear level by using the 
circularly symmetric ansatz. Another extension that is done for the flat spacetime is adding the 
gravitational Chern-Simons term to the general quadratic theory. With this addition the NMG 
theory is investigated and in this limit it is found that the oscillators decouple with different masses. 
Therefore, this model is parity violating theory as expected. 

In dS spacetime case, the general action is written in terms of three gauge-invariant functions. 
These functions are constructed from the derivatives of the components of metric perturbation. 
The fields are decoupled by two different ways. The first way to decouple fields is to go to the 
nonrelativistic limit by dropping out the two dimensional Laplacians in the action. The second 
decoupling way is to get the field equations in the relativistic form. 

These gauge-invariant actions that are constructed in this chapter may be useful when nonlin- 
earities and interactions are introduced to the theory. Apart from this usefulness there are other 
interesting points about the theories that we discussed in this chapter. Tuned values of the param- 
eters give rise to uncommon phenomena, for instance partial masslessness or chiral gravity, which 
especially arise in (anti)-de Sitter spacetimes. These matters are open to study. 

Part IV 

Conclusion 

In this dissertation, the most general quadratic curvature massive gravity theory is analyzed to find 
a specific unitary theory. For this aim the one-particle scattering amplitude is calculated. Also, to 
understand the found unitary theory, that is NMG, the most quadratic curvature theory is written 
in canonical form in three dimensions for both flat and constant curvature spacetimes. 

In the second chapter we have studied the D dimensional quadratic gravity theory augmented 
with a Pauli-Fierz mass term. First the linear theory without sources is studied for both massive 
and massless limits. For the massive case, a wave equation is found with two massive excitations 
in flat spacetime. In curved spacetime the trace of the metric perturbation becomes a dynamical 
scalar field. When the dynamical part is dropped by fixing the coupling constants one solution of 
this equation is the partially massless point which comes out only in curved backgrounds. For the 
massless limit both the linearized Einstein tensor and the linearized Ricci scalar are background 
diffeomorphism invariant. Also, the scalar curvature must be zero in flat spacetime, but need not 
be zero in curved background when the dynamical part of the equation is eliminated by fixing the 
constants. 

After the analysis in the linearized level, we have moved on to compute the one-particle scat- 
tering amplitude between two covariantly conserved sources. To find the exchange amplitude the 
perturbation part must be expressed in terms of the energy-momentum tensor. However, the com- 
ponents of the perturbation part of the metric are not independent. Hence, it is decomposed into 
its parts to obtain the independent components. After computing the necessary elements, the scat- 
tering amplitude is calculated. From this general amplitude equation, the poles and residues can be 
calculated to have an idea about the particle spectrum of the most general theory. Moreover, the 
Newtonian potentials can be computed from this amplitude equation. Nevertheless, the equations 
are very complicated for the most general case. The calculations are done only for some interesting 
limits. First the massive theory is considered in flat spacetime. For this case it is seen that the 
theory suffers from ghosts, and to get rid of the ghost term the coupling constant of the square of 
Ricci tensor must be taken to zero. Also, the Newtonian potentials are calculated for static point 
sources in four and three dimensions. From the four dimensional Newtonian potential, the vDVZ 
discontinuity is obtained. The massive gravity theory does not fit in with the pure Einstein theory. 
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But in three dimensions it works well with the pure Einstein theory for small separations of the 
sources since the massless limit does not exist. The other interesting limit is the massless theory 
in flat spacetime. In this case first mass is set to zero and then the flat spacetime limit is taken. 
For this case the poles, corresponding residues and Newtonian potential energies are calculated. 
When the theory is constrained to ghost and tachyon freedom, from the poles, residues and New- 
tonian potentials, it is seen that the dimension must be set to three, the coupling constant of the 
Einstein-Hilbert term must taken to be negative and there must be a relation between the coupling 
constants of the higher curvature terms as 8a + 3/3 = 0. This theory is known as NMG. Therefore, 
this theory is obtained from a different perspective than that of |35| [36] . In three dimensions the 
Newtonian potential energy again has the same potential energy as the Einstein gravity at the 
NMG limit. The Newtonian potential energy in four dimensions has a ghost term and to get rid of 
this term /3 — ^ limit must be taken. Also, for this limit the general Newtonian potential energy 
is calculated for D dimensions. From this equation it can be seen that for all dimensions there is 
a ghost term. 

In the third chapter, the three dimensional quadratic curvature theory is studied in more detail. 
To understand how the NMG theory is singled out among other three dimensional theories, the 
most general action is written in canonical form in both flat and dS spacetimes. First the flat 
spacetime case is studied. To get the canonical form, the metric perturbation is decomposed in 
terms of six scalar functions. From these functions three gauge-invariant functions are constructed 
by the help of Bianchi identity. The components of linearized Einstein tensor is written in terms 
of these gauge- invariant combinations. With the help of these components and the linearized 
Ricci scalar, the most general quadratic action is written in the canonical form with these three 
gauge-invariant fields, one of which is automatically decoupled from the others. The decoupling 
of the other two fields depends on whether 2a -|- /3 is zero or not. These two cases are discussed 
separately. When this constant is set to a generic value apart from zero, one of the fields that 
has no dynamics in the action can be written in terms of the other dynamical field. After this 
manipulation, a decoupled action is computed. This action is analyzed for some special points 
and it is found that for NMG case, the higher-derivative term disappears and we are left with a 
parity-invariant theory. Except for the NMG case, this theory describes a ghost like excitation. 
When we set 2a -|- /3 = at the action level and define a new field, a decoupled action can be 
written which has a ghost excitation. Then the discussion is extended by adding static sources 
and spinning masses to the action separately. For the static case an interaction part comes to the 
theory. The Newtonian potential energy of this interaction part becomes attractive for negative 
K. By adding spinning masses, the spin-spin interaction part is also found. Unfortunately, this 
interaction can be repulsive or attractive depending on the signs of the spins. Moreover, some of the 
found results are also obtained from the nonlinear theory for the NMG limit. The last part for the 
flat spacetime discussion is to add the gravitational Chern-Simons term to the general quadratic 
action. For this case the masses of the excitations are found. Since, these masses are different, this 
is a parity- violating theory. Also, for the topologically massive gravity limit, the expected result 
is found. 

In the second part of the third chapter, the discussion is extended to curved backgrounds, 
namely to dS backgrounds. In this part the calculations are repeated for this case and again three 
gauge invariant functions are found. With the help of these functions the generic action is written 
with one decoupled field and two coupled fields. For the decoupled field the mass is computed, and 
for the coupled fields the mass is obtained from both the nonrelativistic limit and the equations of 
motion. 

In the Appendix A, we give some details of the computation which can be helpful to follow the 
discussion. 
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Appendix A: Metric Decomposition and The Higher Curvature Action in dS Background 

1. Introduction 

In this chapter some of the calculations are given to follow the bulk of the dissertation. Here 
the gauge invariant action is calculated in dS spacetime. For this purpose, first the gauge invari- 
ance of the metric perturbation is written, then from the components of the Einstein tensor the 
gauge-invariant functions are defined. Then the action is written in terms of these gauge-invariant 
functions. In this part all the raising and lowering operations are made by the background metric 
and all covariant derivatives are defined with the background metric. 

The metric is taken as 

9fiu = Qfiu + hf,u, (Al) 

where 5^1, = ^rj^iy is the background metric and ry^j^ is the flat spacetime metric. The perturbation 
part of the metric /i^jy is decomposed as 



hi 



hoi = {-eijdjTj + d,NL) , /ioo = N. (A2) 



Here Cij is the tensor density with the convention €12 = 1. 
The trace of the perturbation metric is 



h = g^^h^u = ^J?^'' V = -9 (-^00 + ha) 



a2 



{-N + 4> + x)- (A3) 

Some useful identities follows: 

dihoi = di {-CikdkV + diNi) = dfNi, 
didjhij = didj (Sij + didj^ cf) - didjX + {(-ikdkdj + ejkdkd^ i , 

= di [{di -di)(l) + dix - Ukdki] = dlx, (A4) 
ha = [bii did^ (j) - didix + (eikdkdi + eikdkd^ ^ = + 

Later the six free functions are redefined as ■0 = j^ip^ to get simple form of the results. With 



this redefinition, the metric (Al) takes the form 



^2 

= ^ (^/^!^ + hf,u) ■ (A5) 

In this appendix only the calculations of the dS part is given but the flat spacetime part is given 
at the needed steps by taking the limits >-oo, ^/t— >■!. 

2. The Gauge-Invariance of Metric Decomposition: 

The gauge transformation is S^^hf^^ = ^fiCu + ^uCfi where 

Cm = (Co, -(ijdjC + diK) (A6) 
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and eij is the Levi-Civita tensor density with the convention ei2 = 1- To calculate the gauge- 



invariance of the metric components the connection coefficients are needed. For ( Al ) the coefficients 
can be calculated by use of 



(A7) 



where for notational simplicity we omit the bar sign on F^^ and in what follows we shall also omit 



the bar on the covariant derivatives. There are only three non-trivial components of (A7) 



00 



t 



(A8) 



All other components of the connection coefficients are zero. For the full tensor the gauge invariance 
can be written as 



after decomposing the covariant derivative. 



h',.„ - h„,u = d^Cu + a,C/. - 2F^,Ca, 



(A9) 



a. The /iqo Component- 



Let start with the /iqo component. Using (A9); 



^00 ~ ^00 — 29oCo — 2FooCo — 2FooCi 
N'-N = 2Uo + ]Co] , 



where in the second line, (|A8|) is used. Taking /i^^^ ^^/j-i^ 



^(N = 2^-, { Co + ko 



(AlO) 



(All) 



The hiQ Component: 



The hio component can be transformed again using (A9) as 

Ko - hio = doCi + diCo - 2F^oiO - 2FO,Co 

2 

-eijdjTj' + diN'^ + eijdjT] - diNi = do {-eijdjC + d^K) + diCo + -Sij {-ekjdkC + djK) 
-ei.dj (??' {N'l - Nl) = -dje,, U - ^c) + di (k + (o + , (A12) 



Equating both sides with respect to their coefficients, we found 



where we have used (A2) and (A8) in the left and right hand sides of the second line, respectively. 

(A13) 



6^Nl = ( k + Co + , S^V=U-lc], 



and for the redefinition of the metric perturbation, (A13) take the following form 

k^L = ^ ( + Co + j , -^c^ = ^ ( C - -C 



(A14) 
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The hij Component: 



From (A9) the last term is can be written as 

= i^ij + didj^ (j)' - BidjX' + {(-ikdkdj + ejkdkd^ C' 
- {^ij + didj^ (j) + didjx - [(^ikdkdj + ejkdkd^ C 
=dj {-eikdkC + diK) + di {-ejkdkC + djn) + j6ijCo, 



where we have again used ( A2 ),( A6) and (A8 ). The (A15 ) is multiphed with dj in order 
(j) field. Note that djdidj = —di and didi = — 1. Then 

2 

diX' - (^ikdk^' - diX + eikdki = -eikdjdkC + 2didjK + -diCo, 
and equating both sides with respect to their coefficients, we get 



di{x'-x) = 2d]K + -Co, 



5a = 2[ d^K + -Co ) , 



and 



-eikdk iC' - C) = - ^ikdjdkC 
=d]C- 



For the redefined metric, (A17) and (A18) becomes 



1 



+ tCo 



and 



We can eliminate ^ field by multiplying (A15) with and using (A17) we get 



and for the redefined metric perturbation (A21) reads 



For the flat spacetime case the non-gauge-invariant functions become 



0, b^i = d% s^x = 2d]K, 



5^Nl = k + Co, 6^r] = C, S^N = 2Co. 



(A15) 
to eliminate 

(A16) 



(A17) 



(A18) 



(A19) 



(A20) 



(A21) 



(A22) 



(A23) 

From (A23) it can be seen that apart from (j) function all other functions are not gauge-invariant. 
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In this part we see that the functions are not gauge-invariant. Therefore, the next section is 
devoted to find gauge invariant functions. For this purpose, we need to find the components of the 
hnear Einstein tensor in terms of the free functions since it is known that hnear Einstein tensor 
must satisfy the Bianchi identity. Therefore, its components must be gauge invariant. 

In the next section we first find the components of Ricci tensor and scalar in terms of these 
six free functions. Then using these components, the components of the hnear Einstein tensor are 
written. Then we define the gauge invariant functions and also check their invariance. Moreover, 
we see that Bianchi identity give us a relation between these gauge invariant functions. 

3. Gauge-Invariant Combinations: 

The Bianchi identity gives a clue to finding these gauge invariant combinations, constructed 
from the six free functions, since ^ = 0. Therefore, the components of the linear Einstein 
tensor must be composed of such combinations that are invariant under a gauge transformation. 
The components of the linear Einstein tensor are found one by one. 

Before calculating the components of Einstein tensor the components of the linear Ricci tensor 
and the Scalar curvature must be calculated. 

a. Ricci Tensor: 



The Ricci tensor is 



1 



Extracting the covariant derivatives, the Ricci tensor becomes 



(A24) 



and 



R 



1 

2 



+ 



~ 9 ^ dp (j^ ^u^Xcr^ ~ ^ ppJ^au^Xcr ~ ^ pu^ na^>^cr ~ ^ pa^ fiu^Xc 



(A25) 



(A26) 
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after collecting terms in parenthesis we obtain, 

—2 dp ^r^;//iAcr^ ~ ^ pfj^ au^y^tj ~ ^ pi/^ pa^>^cr ~ ^ pa^ pu^Xa | 

- [dpd.h - r^p.dxh) . 

Putting = ^ri"P and summing a and p indices we get 

|9o {dph„o + duhpo - d^h^^) - V^^ {daKo + d„hao - doha„) 
~ ^Ou {dphaO + dahpo — dohpa) — Too {dphiya + d^hpa — dahf^iy) 



R 



pv 



t' 

'2F 



t 



2£2 {dp,hui + dyhpi - dihpy) - Tf^ {daKi + d^hai - dihau) 

^iu (.9phai -\~ doihp^i dihpoi^ i^dphyo^ + d^/hpa doihp^i;) 



^« pu^Xij ^ip^au^Xi ^iu^ pa^Xi ^ii^ pu^Xa | 



l[Op 



d^d,h-T^^,dxh] . 



With the help of (|A8|, (|A28j) becomes 
+2 



t ( 2 

-^M!/ ^ 1 (t^M^j^o + duhp^o - dohp,y) + - {dp.Ko + d^h^o - dohp^) 



+ - {dpKo + d^hp^o - dohp„) 



+ ^ <{ 9i {dpKi + d„hpi - dihp^y) - {do,Ki + dyho,i - dihc 



{dphai + dahpi - dihpa) + - {dpKo + dyhp^Q - dohpi,) 



-2 di (r^j,/iAi) - ^fpTaiyhxi - r° r^^/iAi + -'^pu^XO 
\ (dpd.h - Tl.dxh) , 



(A27) 



(A28) 



(A29) 



finally, we end up with 



R 



111/ 



2£2 



do {d^Ko + duhf,o - doh^u) - di {d^Ki + dyh^i - dih^y) 



With this result the components of the Ricci tensor can be calculated. 



h. The Rqq Component: 



Using (A30), i^oo component becomes 



-"■00 



2^2 



{^0 {dohoo + dohoo - c?o/ioo) - di {dohoi + dohoi - dihoo) 
1 



+ - {dohoo + dohoo - doho) 
+ (dahoi + dohai — dihao) + {dohai + dahoi — dihoa) 
-2 do (rgo/iAo) - di (T^Qhxi^ + -TQ^hxo + TiQT^Qhxi + Tf(,T^^hxi 

^ (dodoh - vl^dxh 



rearranging and doing some cancellations yield 



00 — ~2£2 1 ^o^o*' ~ '^dodihoi + di hoo + -dohoo + 2r"Q {dahoi + dohai — dih^o 



-2 do (rgo/iAo) - di {VQ^hxi^ + -T^Qhxo + 2TfQT^Qhxi 
\ (dodoh - r^o^A/i) . 



Summing the repeated indices and using ( A8 ) gives 

f2 



-n-00 



(dohoo - 2dodihoi + dfhoo^ - ^ (39o/ioo - ^dohu - 2dih 



Oi) 



^ (dodo + ^do ] h. 
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Putting (|A2l) in (|A33|) yields, 



00 



~2P 

2 If 1 



(-N + ct> + x) 



2£2 



(at + a^AT - 2a2iVi) - ^ (sat - 20 - 2x - 252 Ai) + ^ (<A + x) 



2 

f2 



+ 



(iV-0-x) + ^(A-0-x) + ^(A-0-x) 



2£2 



and after the cancellations the final answer of Rqq component becomes 



00 



2£2 



[dfN - 2dfNL + + x) + ^ (2A - 3(/< - 3x + 2d^NL) + ^ A. 



(A34) 



(A35) 



Also, taking the tjl ^ \ and £ — t- 00 limits the flat spacetime case of this component can be got 

«oo = -\ (plN - 2djNL + + x) ■ (A36) 



c. The Rq^ Component: 



From (|A30l) term can be written as 



R 



2£2 



{do {dohjo + djhoo - dohoj) - di {dohji + djhoi - dihoj) 



t 



{dohjo + djhoo - dohoj) 



+ Tio {dahji + djhai — dihaj) + r^- {dohai + dahoi — dihoa) 

—2 do (Toj^Xo) ~ (Tojhxi^ + -Foj/iao + ^To^aj^M + ^ij^Oahxi 



(A37) 



rearranging and doing some cancellations we get 

f2 



L * f 2 1 

Roj = 1 dodjhoo — dodihij — djdihoi + hoj + -djhoo 

+ TiO {dahji + ^j/loj — dihaj) + rg- {dohai + Sa/loi " dihoa) 

-2 do (Tojf^xo) - di (Voj^Xi) + -rpj/iAo + Tforaj^Ai + ^^^^0^/1^1 



(A38) 
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d. The R^f, Component: 



The last term is i?^^ and it can be written from (A30) as 



Rjk = {9o {djhko + dkhjo - dohjk) - di {djhki + dkhji - dihjk) 
+ ^ {djhko + dkhjo - dohjk) 

+ r°- {dahki + dkhai - diKk) + r"fc {djKi + dahji - dih. 



(A45) 



-2 do (Vjkhxo^ - di (Vjkhxij + -Tjkhxo + ^ij^ak^M + ^ik^, 
Ud.dkh-T^kdxh 



and using (A8) yields 



R^h = (2dodjhko - 2djdihki - dohjk + dfhjk 



- ^ C^djhko - Sdohjk + 2r]jkdohoo - 2r]jkdihoi) 



+ 'fi2^jk - 7^ djdkh + -rjjkdoh . 



Putting (A2) in (A46) gives us 



(A46) 



R 



jk 



2P 



d. 



2ekidi (rj+]r]]+dk(2NL-N+^NL 



di + ^^do + ^ 



4 

t2 



(djk + (/> - djdkX + (ejididk + ekAdj^ ^ 



+ 



5jkdf4> - {ejididk - ekAdj) ^ 



1 , 

+ -Ojk 



2{^N-dtNL) + [-N + ^ + X) + -{-N + <P + x) 



(A47) 



where {cjididk — ^kididj) ^ should be eliminated due to the fact that j and k are symmetric. Then 



the final form for R^^ becomes 



R 



■jk 



2P 



1 



2ekidi {rj+-rj] +dk{2NL-N+-NL 



-{di + ho + ^. 



+5jkdi (/) + j6jk 



(6jk + djdk'j - djdkX + [(^jididk + ekAdj^ C (A48) 
(iV - dfNL) + (-N + ^ + x)+j{-N + (t> + x) 
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The flat spacetime result of this component is 



^jk 



2ekidi {r]+^i]]+dk{2NL-N + -Nl 



dl + "-do + 



-6jkdf(j) + -6jk 



1 
t 



{djk + djdkj (f) - djdkX + {(-jididk + ekididj^ C 
2 (iV - dfNL) + (-iV + + x) + ^ (-iV + + X) } , (A49) 



where the i/£ — t- 1 and £ — t- oo Hmits are used. 

In order to calculate the components of Einstein tensor the Ricci scalar is also needed. The 
next section is devoted to the calculation of Ricci scalar. 



4. Ricci Scalar 

Now, let us work on the Rl term: It is defined as 
Rl = ^aV^h^" -ah- 2A/i, 



(A50) 



where in the second line the upper indices are lowered. In (A50) there are two terms that must be 
written in terms of the free functions. These are 



(A51) 



and 



VpVah^u = dp iV^h^u) - (V„ V) - iy.Ku) - T^py (Vahpa) , 

= dp (j^ah/Mu — ^ (jphxy — T ^^hpxj — T (j^ahpu — r^j^/iAi^ ~ ^au^ 



(A52) 



After the application of contractions, using (A8) and doing some simple manipulations these terms 
become 



nv 



£2 



^2 



dodoh 


- didih + 


\{doh 


dodoh 


- didih - 


Idoh 


{dodoh 


- didih) - 


^doh 



(A53) 
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and 

£4 



74 c^o (Sq/ioo — Too^Ao — ^oohoxj — Tqo (dahoo — r^o^Ao — r^o^oA^ 



00 



Too^Aa 



Toq^OA 



-r 



00 



do (dihoi — T^q/iaj — T^ihoxj — rgj (dahoi — T^g^Ai — T^j/ioA 

[diha 



di (dohio — TQ,-hxo — T, 



^ h 
00' HX 



^oi (d'iho, 



^iohxa — ^iahoX 



^ai^XO 



r^o^jA 



~^fi (pohaO — Toa/iAO " ^QO^aX^ — ^to (pohia — ^Oi^Xa — ^Oa^^iX 



di (d. 



"pA pA 

ji jj 



^ai^Xj ^aj^iX^ 



Summing the repeated indices in (A54) and again using (A8) yields 

f3 



{dodohi 



00 - 2dodihoi + didjhij) + (Sq/ioo 



dnh 



'or Hi) 



Putting ( |A53[ ) and ( |A55[ ) into (A50) and using K = l/r gives us 

+4 



Rl 



(5o9o/ioo - 2dodihoi + didjhij) + {dohoo - doha) 



+ 



{dodo - didi) 



h. 



(A54) 



(A55) 



(A56) 



Putting (A2) in (A50) yields 



Rl 



^N-2dfNL + dfx) + ^^{N 



X 



f t 
j2 i^odo - didi) - j^do 



£2 



{-N + (t> + x) 



(N-2dlNL + djx) + y,{N-^ 



X 



4^3 

{dodo - didi) + --do + 



2f 



£4 £4 



do 



2t^ 2r 



{-N + ^ + x), 



j-,dodo + j-,do]N + 



£4 



d 



do\x-^,2djNL 



do(, 



3^3 

^ {dodo - didi) + —do 



2*2 

1^ 



j^d^di 



2t^ 2*2' 
^^0 + ^ 



N + 



2*3 



2t^ 



X ■ 



2dfNL 



+ 



*4 2*3 
^ (ao9o - didi) + — 9o 



2*2 



(A57) 
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and finally the scalar curvature read 



RL='ji{dfN + x-2dfNL + ^-dfcl,) + '^[-N + ^ + x)-'^{-N + cP + x). (A58) 



2t^ 



In the flat spacetime limit ( A58 ) becomes 



Rl = d]N + X - 2dfNL + {di - 9: 



(A59) 



which is the flat spacetime result for the scalar curvature. 



With the results that we found in the above calculations, (A35), (A43), (A48) and (A58), we 



are ready to flnd the components of the Einstein tensor in terms of the free functions of metric 
perturbation. 

5. The Einstein Tensor 

The linear Einstein tensor is deflned as follows 



Ri 



The components of the Einstein tensor is investigated term by term. 



(A60) 



a. The Qkr, Term: 



The 00-component of linear Einstein tensor, ^qq, can be written by using (A35), (A58) and (A2) 
as follows, 



^00 — -^00 ~ 2^00 Rl - 2 A /loo, 



(9f iV - 2dfNL 



+ <P + X + 



2£2 



(2N 



3x + 2dtNL) + -^N 



dfN + x- 2dfNL + {di - 5| 



-—[-N + cP + X 



2t^ 

1^ 



i-N + cP + x) 



'2P 

' f_ 
2P 



(dfN - 2dfNL 



+ <P + X] + i^\2N 



3X + 2d^NL 



dfN + x- 2dfNL + {di - df) + 1 (-AT + + X 



1 



{-N + 4> + x)j, 

(A61) 



and doing cancellations and summations we get 



^00 



(A62) 
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To simplify (A62), the metric perturbation functions are redefined as 



and the time 



derivative transforms as 15 
Einstein tensor becomes 



-2|g-0 + With these redefinitions the 00-component of the 



^00 



t 



-2j^ct^+-^-2^X+^X-2^dlNA--,{-N + 4> + x) 



1 : 1 



+ Y^X-\d^Ni)j-^^[-N + <t^ + x) 



t 



dfNL 



From the last line of (A63) a gauge invariant function can be defined: 

1 



/ 



t 



+ 



then 



iV2 



^00 



+ X--N 



t 



Nl 



2i 



Also, for flat spacetime this term becomes 



(A63) 



(A64) 



(A65) 



^00 



(A66) 



where f = (j) for the flat spacetime background. 



b. The Goj Term: 

The Oj-component of linear Einstein tensor, ^q^- = Rqj — 2A/ioj, can be written by using (A43) 



and (A2) as follows. 



r 

'2P 



dj(j) + ejkdk (i - dfrj 



dj4> + ejkdk (i - dfri 



t 



+ i^[dj{N-24>)-2ejid,i\. 



(A67) 



Again doing the same redefinitions 

2 



^4> and 



-2^4>+ ^4>, (|A67|) becomes 



1 



+ - [dj (N - 2<P) - 2ejidi^] 



From ( A68 ) two gauge- invariant combinations can be defined as 

£ 



P 



t 



t 



t 



(A68) 



(A69) 
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With these definitions, (A69), the Oj-component of the hnear Einstein tensor reads 



21 



For the flat spacetime this component yields 

%j - "2 

with the flat spacetime version of the gauge-invariant functions, that are 

p = 0, a = 



dj(j) + ejkdkcr 
ant functi 



(A70) 

(A71) 
(A72) 



c. The Q^j^ Term. 



The jfc-component of Hnear Einstein tensor, can be written by putting (A48), (A58) and 



( A2 ) into the following equation 



that yields 
Qjk = - 



212 I 1 



Gj"k — Rjk - 2^^ikRL - -p^hjk, 



-2€kidi {ri+jv)+dk (2Nl -N + ^Nl 



(A73) 



t2 



(6jk + djdk'j - djdkX + [(^jididk + eudidj^ i 



-Sjkdi (/) + jSjk 



«9 ""7 A; ) 



2 {^N - dfNL) + [-N + cP + x)+-{-N + ^ + x) 
2t 



(A74) 



^ [dfN + x-2dfNL+{di-df)cl>]+'^[-N + ^ + x)-l{-N + ^ + x) 



and with some cancellations it takes the following form 



2(2 I 1 9j 



-2ekidi {ii + ^v)+dk (2Nl -N+'^Nl 



5o + fo 



djdk4> - djdkX + {(-jididk + ekAdj^j ^ 



\ ^ {dfN + x- 2dfNL) + 2dfNL 



3t 



(A75) 



£2 



After redeflning all the functions as 

j4 f3 ^ + fi4> respectively, (A75) becomes 



for which the flrst second derivatives are 



-^(p+^(f), and 



-2gff, = 2ekidjdi [^v-v]+ djdk (-'^Nl + 2Nl-N 



+ djdk ( ^< 



djdk [^X-xj + {ejididk + ekAdj^ - | 



(A76) 



+ 5jk { dfN --X + X- 2dfNL + ^dfNL + - 



2 



1 V, 
t 
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Since the j and k indices are symmetric, the first term can be manipulated such that 
2ekidjdi [jT] -f]j= {ekidjdi + ejidkdi) Qr? - ?? 



1 



(A77) 



where in the second hne the unit derivatives are introduced. After this manipulation (|A76|) becomes 

->L 



-2g^k = - (ekid.di + ejidkdi) (^^dfii - dfri + C 
-djdk(^-j4> + 4'+^N-^N 

+ {Sjk + djdk) (dfN -jx + x- 2dfNL + jdfNL + ^N-^N 



(A78) 



where the terms ^dkdjN and jdkdjN are added and subtracted. The second parenthesis in the 
first term is & and the second term is p. The last term can be defined as 



^ (dfN --^x + X- 2dfNL + -^dfNL + ^N--^N], 



(A79) 



which is also gauge-invariant function. Finally, the j/c-component of the linear Einstein tensor 
becomes 



_ t 
'^'^ ~~2£ 

The flat spacetime version of this component is 



^jfe = {- {^kidjdi + ejidkdi'j a - djdkP + (Sjk + l} 



Qjk = -\{- {ekidjdi + ejidkdi^ a - djdk4> + {^jk + djdk) q\ , 



where the function q reads 



q = dfN-^X + X-2dfNL. 



(A80) 



(A81) 



(A82) 



6. Gauge-Invariance of the Defined Functions 

The functions can be checked whether they are gauge-invariant or not. Let us start with, 

S^a = J [6ci - V^6^r^) , (A83) 



and using (|A11 A14 A19l A20 A22|) and their derivatives (|A83|) becomes 

0, 



i 2t, 



2t, 



(A84) 



where we have used (5^^ = ^djC + ^djC- Again using (All 
reads 



A14 



A19 



A20 



A22 ) the p function 



5(^p 



t 

1(2 





1 






2t. 






-4 



1 



(A85) 
(A86) 
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where we have again used the derivative of fields, S^cj) = + ^Co- The other two functions can 
be handled in the same way. 



+ 



1 

tV2 



+ - -k^ - -k^L 



(A87) 



tV2 

2t 



2t ■ 



t 



^2^^ ' £2' 

2t 



1 



Co + T^Co + 4^ V^K + tCo + V"/i + -Co - T^Co 



t 



f2 



1 



1 



t 



Co + -Co^ 



£2<0 ^2 



+ Co + ) =0, 



(A88) 



where we have used = 4^ (^j + iCo) + '^ja {d'jk + jCo — ^^Co) • The q function is 



t 



ki 



2*2 
4 



Co + ^Co 



1 



V^K + ^Co 



4t 



4t 



2t. 



t 



1 • 1 



V^K + tCo ) + ^ ( V^K + ^Co 



1 



t 



t2 



2t 

Co 



v^K + -Co - ^Co 



t2' 



2t' 



+ ^ ( v^/^ + Ko - T^Co ) + ^ ( V^Ai + -Co - ^Co - -^Co + ^Co 



1 r- 1 • 1 • 



t 



4 



V" ( K + Co + 

2 

-/ 
t 

2t 



2t^ 



2 2 

( K + Co + -K - 



+ 7«V^ K + Co + T^ +^Ko + tCo 



Co + ^ Co 



Co + -Co - ^Co 



0. 



(A89) 



(A90) 



Therefore, these functions are indeed invariant under gauge transformations. 



7. Ricci Scalar in Terms of Gauge Invariant functions 



The Ricci scalar can also be written in terms of the gauge-invariant functions. We first write 

). With these functions the Ricci scalar takes the form 



(A58) with the modified functions. 



£2 



d? 



0/1 O f-1 o 



(A91) 



Using the definitions of the gauge-invariant functions, (A64), (A69) and (A79), the Ricci scalar 
can be written as 



i?L = ^(<?-VV + p). 



(A92) 



For the flat spacetime case 



(A93) 



and using □ = —do + V'^ the Ricci scalar becomes 



RL = q- □</>• 



(A94) 
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8. The Bianchi Identity 



The Bianchi identity can also be written in terms of the gauge-invariant functions. First the 
identity is decomposed into its components by doing the summation in the repeated indices: 



— 9' y^cr^^u 

9 







(A95) 



Xil ) 



£2 



Of 



rO /^L -pU /^Ij "p/i: 







(A96) 



There is two equations depending on the u index. For = 0, (A96) yields 

+2 



+ 



£2 



L -pO /^L pO /^L -pk pL 

^ 00*^00 ^ 00*^00 ^ nny, 



OO^feO 



L pO /^L pO /^L -r^k /->L 
iO ~ -"^ jif 00 ~ jOfOi ~ J- jO^fci ) ' 



(A97) 



and using (A8) 



+ 







£2 



50^00 - diG, 



(A98) 



Using the gauge-invariant form of the components of the Einstein tensor and their derivatives, that 
are 



and 



^00 



t 

U 
t 



21 ■' 21 



t 

'2I 



dfp 



21 



(A99) 
(AlOO) 



21 
t 

'2I 



{p + q)- 



(AlOl) 



With these equations, (|A99|), (|A100|) and (|A101|), the Bianchi identity yields 

t 







1 



* -v2i , * V72. 1 



21 ' 21 ' 21 21 
[-V^f-tV^f + tV^p + p + q) , 



(A102) 
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and finally 



I ^f + f-p] -p-q = 0. 



Also, u can be taken as u = n. For this choice (A96) reads 

+2 



= --^ 



"pU nL -pU -pfc cL 

^ OO^^On On»00 flrv.f fci 



On»A;0 



nO 



1 







^oQL-\<^L-^^Ql^ 



t 



Again using (AlOO) and the derivative of (|A80|), that is 

diGi 



-•L 



t 



21 



{-(^nidia + dnP) , 



(A104I yields 



(A103) 



(A104) 



(A105) 



[dnP + enkdkcr] + [dnP + enkdko] + ^ [-enididia + dnp] , 



2e 



21 



1 



and 



= -dip - eikdko- + - [dip + eikdkcr] - e„/9/o- + dnP, 
= -2e„fc9fccr + ^ [dnP + enkdkcr] , 



= Cnkdk [ -2a + ^cr j + ^dnP- 



(A106) 



(A107) 



Using (A103), (A92) becomes 



Rl = -^[q-V'f+p 



(A108) 



Up to now, we write the components of linearized Ricci tensor and linear Einstein tensor in 
terms of gauge invariant functions. Moreover, the linear Ricci scalar is also written in a gauge 



invariant form. Now we can write (147) in terms of gauge invariant functions by using these 



identities. The next section is devoted to find a gauge invariant action. 



9. Metric Decomposition and Gauge-Invariant Form Of The Higher Derivative Action 



Linearized action is 



(A109) 
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where background metric is 



-df + dx^ + dy' 



(Alio) 



Note that the term in the parenthesis is the field equations of the quadratic curvature theory. The 
constant a in front of Gj^'^ is j27| 



a = - + UAa + 2A/3. 

K 



(Alll) 



Linearized form of Einstein and Ricci tensors, and Ricci scalar are given in (AGO, A24 A50) with 
the definition □ = V^V^ = ^^'"V^V,,. 



The Einstein-Hilbert Part 



Let us work on the Einstein-Hilbert part first: 



— — - / d-^T——r)P'^n'^'^h 



t 



2 I d'x-vP^r^-^h,^gl:,. 



(A112) 



Writing Einstein-Hilbert part in terms of tensor components yields 



(A113) 



Let us calculate term by term this action: 



h. Decomposition of Hq^Q^q term 



Let's start with the first term, that is /loo^ob where 



00 

?ob = -^00 ~ 2^ooRl - 2A/ioD. 



(A114) 



Using (A62) and (A2), (A115) becomes 



hooGno = ^ 



2£2' 



d!cl>-^{^ + X-2dfNL 



j^{-N + (t> + x) 



(A115) 



For the flat spacetime case (A115I becomes 



^00^00 



:N&tct>. 



(A116) 
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Decomposition of h^jQ^, term 



Let us move to second term in Einstein-Hilbert action which is — Using (A67) and 
(A2| hQjQ^j becomes 



-2h 



2£2 



<P--[N-2(t>) 



(djhoj) 



_ _ (AT _ 20) 



{ejkdkhoj) 



{i-dl 



(A117) 



where the integral sign is suppressed and at the last line the boundary terms dropped. Also note 
that, 



and 



Cjkdkhoj = €jkdk {-ejidirj + OjNl) = -ejiCjkdkdir] + CjkdkdjNL, 
= -SikdkdiT] = -dfr]. 



Then (A117) yields, 



dfr] 



<^ _ _ (iV - 20) 



In t/t — ^ 1 and £ — ^ oo limits (A120) becomes 

-2ho,Ql^, = (i- dfr]) dfr] - <P [d]NL 
which is the result for the flat spacetime case. 



(A118) 



(A119) 



{i - dfv) dff] - 4> {d]NL) + J [2Cdfr] + {N - 20) (sJa^l)] | . (A120) 



(A121) 



d. Decomposition of hjkGfk t^rm 



Using (A75) and the first equation in (A2), hj^Q^'f. takes the following form 



hjkQ 



'jk 



hikdi 



-2ekidi {f]+^^r]\+dk (2Nl -N+'^Nl 



2£2 j ] -^i^^J 

-hjk ( dl + jdo ) djdkcj) - djdkX + {(^jAdk + ekididj) ^ 



^hjkSjk \^{dfN + x- 2dfNL) -j^{n + 2dfNL) + 



3t 



X 



(A122) 
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Note that djhjk = dkX ~ ^kidiS, and djkhjk = (j) + X and again the integral sign is suppressed. 
Taking the space derivative in front of the metric perturbation in the first term with dropping the 
boundary terms the equation becomes 



-2ekidi ( 77 + 7^) + dk (2Nl - N +-Nl 



tii) 



-hjk + jdo^ djdk(t> - djdkX + {ejididk + ekAdj) i 

-\{^ + x)^Y^{djN + x- 28^ Nl) + 2dfNL) + • (A123) 

The multiphcation of the metric perturbation with the middle term is 
[Sjk + dkdj^ 4> - BkdjX + (ekndndj + ejndndk) i 

'djdkcj)'^'^ - djdkx^^ + (e.Adk + ekididj) c 



(A124) 



where denotes time derivative and i = 1,2. Putting this equation with suitable derivative signs 



and doing the cancellations (A123) turns to be 



hjkQjk 



2P 



24>dfNL + 2Cdff] + ct>dfN + ct>x + x4>- 



In the flat spacetime limit (A125) becomes 
hjkQjk 



7) + (l)dfN + (t>x + 



2CC , 



32 a: 



and with integration by parts on time differentiations (A126) yields 

1 



hjkQj'k 



^fNL + idfrj + el - ^x 



(A125) 



(A126) 



(A127) 



e. Decomposed form of Einstein- Hilhert Action 



To get the final form of Einstein-Hilbert term, all the found results (A115), (A120), (A125) are 



put in (A113) to get 



t 



(fx-N 



2P 



d}4>-^^(4>^X-'2dlN:)-j^{-N-,4>^x) 



S - ^'^) - <P {Q^Nl) + - [2Cdfri + {N- 20) [d]NL 



dr'x 



2P ) 



] { (-2(t>djNL + 2Cdfr, + ^dfN + cl)x + X^- 2^ 
+j {2^dfr] - - 2(l)dfNL -4>N-xN + 3</>X + 3x0 



(A128) 
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To see the terms that may cancel or sum after integration by parts operation, (A128) is written in 
the following form 



+ l^dfrjf -2idfv + 2Cd!r,-2d 
+^ (3<Px + 3x<A + 2cl)d]NL - 2idlri - 
[n(^^ + x)-N{^ + x)- 4NdfNL 



+- 



+-Ni-N + cP + x) 



After doing integration by parts in the terms x4' N ((j) + x): (A129) becomes 



2^3 



6 6 

<PX + ^[X + -X+ ^X]+ '^^d'fN 



2(PdfNL - "-(pdfNL - 2cl)dfNL 



6 



^2 aV 



+ ^ (s^X - 3X0 - ^X^ + 2cl)d]NL - 2Cdfrj - 6^ 



+ - 



N + x) + N + x) + {(j) + x) - ^NdfNL 



+ -N{-N + ^ + x) 



and the final form of (A130) in terms of metric perturbation functions is 

J Ux + <t>df (n - 2Nl) + {-dfv + if 



J (30X + n[<P + x)- 2NdfNL - 2(t>dfNL - 4^9^ 
_L(_Ar2 + 2iV(0 + x))}, 



which is not gauge invariant yet. In flat spacetime limit (A131) yields 



4>x + <l^df [n - 2Nl) + {-dfv + if 



(A129) 



(A130) 



(A131) 



(A132) 
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/. The Gauge-Invariant form of the Einstein- Hilbert term: 



From (A131) the gauge-invariant form of the Einstein-Hilbert action can be written. For this 



aim, first the functions are redefined as (/) — >■ ^(p. With this redefinition (A131) becomes 
a 



Ie 



t 



+ 



<t>x + '\>dl (iV - 2Ni) + [-dh + i 
1 + AT + ^) _ 2NdlNL + 24>dlNL 



._JV2 J £2 



and doing integration by parts in the last term 



(A133) 



(A134) 



where we have dropped the boundary term in the last equality. Then the Einstein-Hilbert action 
becomes 



1 



t2 



4>x + (iV - 2iVi) + (-dh ^i) - 



(A135) 



where the integral sign and the overall coefficients are suppressed. Using (A69) and adding and 
subtracting the terms f (/'A, 0^, yields 



Ie 



^2 



0X + (30 + x) - 2A52a^l + 2</.a2 Al - 



(A136) 



Again doing integration by parts with respect to time in the following terms terms jX^; -j^'^ ^^l 
and adding the term t0V0 freely, since it is a boundary term, to the action yields 



Ie 



+ - 0x + ^dl^' + 209f A^L - 72 - i0V0 



i2 



2A0 A0 Ax - 2Nd'^NL - 0^ 



(A137) 



With reordering the terms (A137) can be written as follows 

2 



Ie 



a'+p')-t4>[\x--V^NL-'^N+-^4> + ^' 



t 



,1 2 9 2 1 • n 

N[-x--V^Nl-^N+j<P + V^ 



and putting (A64) in (A138) we obtain 



(A138) 



(A139) 
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Finally, using (|A69|) we get 



Ie 



(A140) 



We can also do some manipulations to write (A140) in terms of the Ricci Scalar. When we look 



at the gauge-invariant Bianchi identity, we see that the term tfV^f + fV'^f must be added to 



the action (A140). Since this term is boundary term it can be added freely. Therefore, the action 
becomes 



and using Bianchi identity (|A141|) becomes 

t 



lE=[-l) [{'r'+p')+fiq + p) 



Then, we add and subtract /V^/ in (|A142| and using (^A92h we get 

Ie = - 

and in the formal form we have 



iE = ll\l{<r'+p' + fv'f) + pRL 



The flat spacetime limit of this action is 

1 



2k 



(A141) 



(A142) 



(A143) 



(A144) 



(A145) 



where the flat space version of the gauge-invariant form of the Ricci scalar ( A94 ) and a = ^ are 
used. Then 



Ie 



1 

2k 



2k 



(A146) 



where in the first line integration by parts is done and in the second line the definition of the 
D'Alembertian □ = —Oq + is used. Also, at the flat spacetime limit the coefficient a becomes 



after taking A = in (Alll ). 



10. The 2a + (3 Part of The Higher Derivative Action 

The 2a + P of the action is 



(2a + P) 



I d^x v/^V - V^V^ + 2A5'^^) Rl, (A147) 
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and doing integration by parts and dropping the boundary terms the action becomes 



I 



2a+l3 



(2a + /3) f,3 



J y/^Ri {ah - VVh^^ + 2Ah) , (A148) 



where the term in the parenthesis is —Rl- Therefore, the action becomes, 



(A149) 



since Rl is already gauge-invariant this part of the action is also gauge-invariant. The flat spacetime 
version of this action is 



(2a + /3) [ 



d-'xRi, 



where -v/^ = |^ = 1 is taken. 



(A150) 



11. The /3 Part of the Higher Derivative Action 

Now, let us move to the (5 term of ( A109[ ); 

Ir = -^ Jd^x (°^r - A5^"i?L) , 

J d^x ^ [(□ V) - A/ii?L] , 



2 
2 



(A151) 



where in the second line the D'Alembertian is moved in front of the metric perturbation and the 



boundary terms were dropped. Observe that, the definition of i?^^ , (A24), involves ^h^y, therefore 



it can be written in terms of Qj^y From (A24) 



ah. 



and changing the order of the covariant derivatives yield 



+ GAh^u - 2Ag^„h, 



where we have used 



V"V/,/i^^ = [V", V^] h^^ + V^V^/i^^ 



= 3Ah^^ - g^yAh + Vf.Vh^^, 
and in the second line we have put the following identities 

Rlipua = A {g^ugpa — dfiagup) , Rpu = 2Ag^u, R = 6A. 

Then, 



ah 



pu 



-- -2RI^ + 6A V - 2A5^^/i + + V^V^ V - V^V^/i 

= -2g^^ - g^,^RL + 2Ah^y - 2g^^Ah 



(A152) 



(A153) 



(A154) 
(A155) 



(A156) 
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where in the second hne we have used (A60). Then the first term in (A151) becomes, 



(A157) 



The last three terms in (|A157|) becomes boundary term and drops out by the use of Bianchi identity, 



^uGt'^ = 0. Also, note that 



flU 



9^u Rl 



flU 



OfiuR^ — ^Rl 



3 

2' 



2Ah, 



and using 



Rl ^ {g^^R^ 



'fjLU 



r^'R, 

g'^Rflu 



2A/i, 



g'^'^R^, = Rl + 2Ah, 



(A158) becomes 



g/iuG 



-Rl. 



Putting ( |A161[ ) in ( |A157[ ) yields 



1 

2^ 

The total action for the /? part of (A109) becomes 



-2g^',gr + ^Rl + 2A + ^hRL. 



(fx V5 



(A158) 



(A159) 
(A160) 

(A161) 
(A162) 
(A163) 



where A = is used. The gauge-invariant form of the last term is known since it is the same as 
Einstein-Hilbert part. The middle term is already gauge-invariant. To have a total gauge- invariant 
action the first term must be written in that form. 



a. The Gauge-Invariant form of Ga^Gi^ ' 



Decomposing /3 f d^x [pj^yQi^^ gives us 



P I d'x^gj^.g'r = p I d'x 



t^i/^L 



(ao^o)'-2(ao^.)'+(a,§)' 



(A164) 



From ( |A65[ |A70[ |A80| ), we have 



4^2 



2^2 +g2 +p2 



(A165) 



and using (A165) in the action (A164) we get, 

+3 



- t 

gQliuQL = ^ 



2P 

2a^ + g2 + p2 



(A166) 



where ^/—g = ^ and note that we lowered the upper indices by using the background metric. 
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The Gauge- Invariant Form of the Ip Action: 



Adding (A166), square of (A92) and (A144) the action (A163) becomes 

f3 











-IJ 


d^x 






4 


I"' 







1 



(VV) ' + pV^p + aVV + a^ + ]-q^ +'p'^-' (q- V^f + p 



1 



1 



£3 



pV'^p + aV'^a + a'^ -qp + qV"^ f + pV^ 



(A167) 



In this action there is five gauge-invariant functions. By using the gauge-invariant form of the 
number of functions can be reduced to four. Here, it is preferred to ehminate the q variable by 
using = g (g - VV + P) ^ 9 = g^L + VV - p, 



d^x 



pV^p + aV^a + 



t3 



+ 



d-'x 



2 t 
F [i 



^ Jd'x^^ pV^p + aVV + a^- ^^pRl - pV^ + p^ + ^^RiV^f + (v^/) ' 



(r'+p^ + fV'f) + ^fRL 



+ 



/3 



d^^X 



<T'+p' + fV'f)+^fRL . 



(A168) 



This action is the final result for the /3 part of the total action. Now everything is ready to write 
(A109) in terms of the gauge-invariant functions. 



Before going on, the flat spacetime limit of (A168) can be written as 



d^x 



d^x 



d^x 



0V2(/> + (jVV + cj2 - ^Rl - '(jN^cf) + 4)^ + RlV^cI) + (vV)" 



U<j)) - (PV'^cj) + + (g - U(P) V^(P + ( V 



V2 - dl 



-□0 [-dl + 



□c 



(A169) 



where in the second line we have done integration by parts and used ( A94 ) , in the third line some 
suitable combinations have done and in the last line the definition of the D'Alembertian operator 
,□ = —Sq -\- V^, have been used. Finally, for the flat spacetime limit (A168) becomes 



/3 



d^x [aUa + qUcf)] . 



(A170) 



73 



12. The Total Gauge-Invariant Action 



Summing (A144), (A149) and (A168) the total action can be written as 



I 



+ (2a + /3) ^iii 



+ V + p2 + pv^p + (v V) ' - tpRL + S^lV^/ - pV^f 



(A171) 



This result is the gauge-invariant action for the general quadratic curvature theory in three di- 
mensions. This action can also be simplified by defining (p = V^/ and using (A92). For the first 
term 



using the Bianchi identity (A103I 



(A172) 



(A173) 



and observe that '-^fV^f = do{'jfV^f) - f f f - f f = -ffV^f-'^fV^f 
T /^^/ ~ ~J f^'^f after dropping the boundary term. Then (A173) becomes 



t 



a'^ + p^ — t(fp 



The second term of ( A171 ) can be written as 

t3 



(A174) 



(A175) 



where it has been used Rl = fjV^ ^/ — p^ which comes from using Bianchi identity (A103) in the 
gauge-invariant form of Ricci scalar. The last term of (A171I can be written without taking care 



of the a field by using (A108), 



£3 



pV'^p - pq+ (p + q) V^/ 



pV^p -pq + tV^ {f-p)y^f+ (VV) ' 
pV^p -pq- tV^pV^f - (VV)^ 



(A176) 



where in the last line the integration by parts is used for the term l^-V^/V^/ = (V^/)^- 

Again using the Bianchi identity (A103) qp = tV^ {^f — pjp + p i^"^ f) — ■ Then, (A176) can be 
written as 

pv^p - tpv V - tpN'^p - p (vV) + - tV^pV^f - (v V) ' 

f3 



O O 

pV p — tpip — tpV p — pip + p — t(f\7 p — f 



-pV^p — tpip — pip + p^ — tifV'^p — (p? 



(A177) 
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where in the last hne the integration by parts is used for ^pV'^p = — ^pV^j5. The final form of 
dAml) becomes with ^Krf^, (|AT75| and ^nff\ 



I = ^ I d^x {a + 



1 



21 



+ n / ^ ^ InP P ~ ^P^ -p(p + P - tifV P - (f + la, 



where we have defined 



(A178) 



2/3\ t 2 



(A179) 



In the flat spacetime limit, (A178) becomes by summing (A146), (A150) and (A170) 



+ ^ I d^x 



-(Pq + (2a + f3){q- Oclyf + ^qD^ 

K 

1 



(A180) 



13. Gauge Invariance of The Field Equations 

The field equations of the action can also be written in terms gauge-invariant functions. In 
the field equation only the OQj^j^ is unknown in terms of the gauge-invariant functions. Therefore 
in this part only this term is written in terms of gauge invariant functions. First the covariant 
derivatives are extracted 



= ^^V'^' {{dpdaGf:, - GtdpT^ap - r^M^P^i - QpxdpTt " ^LdpGflx) 

(■pQ o -pQf -pA -pa pA \ 
pa'-'aypu ~ pa'- ap^Xu ~ pa'- au^ pX J 

f pQ o pa pA /^L pet pA /-»L \ 

I pp^o-Vau pp'- aa^Xv pp'- av^aX J 

- (Tp^daGf^a - ^pu^apGxa " '^pu^aaGflx) } i (A181) 
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where we have used g^p = jir]ap in the third hne. Doing the summations in the repeated indices 
( A181 ) becomes 



(JV 



I "pa / "pA I -pA /^L \ I / pQf pA i pct pA 

po- I i a^!='A!/ + -L ai^ypX ) ~^ o-a^Ai/ + ^ pp.'- ai/faA 



+ 



/ -pQ; -pA /^Z/ 1 -pa p^A i^Z/ \ \ 

I / p^a p^A /^L I Y^a pA /^_L pct pA , p^a p^A \ 

T" l^"-"^ OAtJ- Oa^Ai/ "T J- i/i-"^ icffAi/ ~ Op'- Oi/faA "r ip'- iv^aX J 

I / p^Q pA /^L I pa p^A /^Z/ p^a p^A /^L i pa pA /^_L \ 1 



(A182) 



and using Christoffel connections (A8), (A182) yields 

f2 



^pu p, 



7L 

^ pu 



3 



'^Qpu^o{-^+Qt^^T^p + Gjlx^^Ti 



pv 



pv 



I ~ nLi I / nh I p^a p^A nh i p^a p^A /^L i opa p^A /^L 



(A183) 



Doing the essential cancellations (A183) becomes 



-(2Tld,Gl + 2Tid,Glix) 

^^fp^^aGxy + ^tv^laGpX + ^rf^r.^^,^^;, | . (A184) 

Again doing summations in the repeated indices and dropping the derivatives of the Christoffel 



connections, diTi^^ = 0, (A184) yields 



□a 



pu 



-dl - ho + df )Gf:,-2 {rldiGi + Ti^diG^ + rld^Gfio + rlAGf^, 



(A185) 



By using (A185) we can calculate the components term by term. 
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The OGnn Term 



Setting = u = in (A185) and using (A8) we have 



□a. 



00 



t 



^ pO r>L ^T* nL ^ pO nL ^T* pL 
+2riori0^ob + 2r°or|o^O^- + ^r^o-'^iO^io + ^r^o-'^iO^iA;} ' 



(A186) 



and after simple manipulations (A186) becomes 



°^oo = ^ 



£2 



-5o^ 



5o+5ne, 



^00 



£2 



-5o^ 



t 



t2 



t2- 



(A187) 



Putting (A65) and its derivatives, the derivative of (A70) and the trace of (A80), which are 



^O^o^O 



-^VV-^VV, «o 



v2/_lv2/, 



'2^ 



L 



t 

'2I 



into (A187) yields 



(A188) 



°^oo 



+ 



2F 



+ 



2F 



t2 



+ 



2P 



(p + ,)-vV 



and from Bianchi identity (A103) the last term is -^Rl and the final result becomes 



t 



2£3 I t 



t2 



(A189) 



(A190) 



For the flat spacetime case (A190) yields 



°^o'i) = ofv'/-V2v2/ 



b. The agj^^ Term 

Let us continue with D^^q component. Inserting ix = k and = 0, the equat: 



;ion 



-dl - ho + Qto - 2 {r^d^g^, + T{,d.Gtj) 



-dl - ho + df^ Gto + \ {duQ^o + 



^2 



'ki I + To^Ok 



Using (A70) and its derivatives, which are 



- ^ {dkP + eknOncr) - ^ {dkP + ekndncr) , 



the derivative of (A80), that is 



dig"- 



ik 



2£ 



(dkP + €kndn<j) , 



and (|A8| into ( |A192[ ) yields 

ff2 



l~ 1 t 

■^^Oko = J {dkP + ekndnCr) + — {dkP + ekndnCr) 
3 1 3 t 

+ - ^ {dkP + ekndnCr) + ^^^"^ ^kndnCr) 

- ^ {dkdfp + ekndndfa) 
2ft t 

5 t 

- ^Y{dkP + ekndnU) , 



and 



agto = ^dk (p + ^p- V^p -Ip- tVV 



t 



2P 



(^kndn \(T + -a - V'fJ - -^a 
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The flat spacetime version of (A196) is 



= \dk {p - V^p) + ]^ekndn (<y - VV) . 
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The OQf^ Term 



The last component can be calculated in the same way. The indices are renamed in (A185) as 



11 = 171 and V = n. After using (A8) 

e 



= ( -di - -do + Of] + - [dmOL + dnOL) 



^ L ^ L ^ L 



(A198) 



After calculating the derivatives and doing suitable cancellations the final answer comes out as 



2^ 



z' 5 1 2 

/' 5 1 2 



(A199) 



In the flat spacetime limit (A199) becomes 



2ag^^ = + d„,dn) [q - V^q) - dmdn (P - V^p 

- {(-mkdkdn + enkdkdm) (^'c^' " V^(T 



(A200) 



With (A190), (A196), (A199) and the results that are computed for Qf^^ and Rl are enough to 
write the gauge-invariant form of the equations of motion. 
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